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ABSTRACT 


The Shuttle Electric Power System Analysis Computer Program (SEPS) was 
developed by TRW under JSC/TRW Task 542 for the Consumables Analysis Section 
of the Mission Planning and Analysis Division. The SEPS program has two 
major uses; first, to perform detailed load analysis including predicting 
energy demands and consumables requirements when the shuttle electric power 
system is operated and perturbed in accordance with premission flight plans; 
and second, to perform parametric and special case studies on the Shuttle 
electric power system. As an additional feature, the SEPS program can be 
and has been used to analyze the ASTP Apollo electric power system. No 
program changes are required to use the SEPS program for analysis of the 
ASTP Apollo electric power system. 

The SEPS Computer Program is written in FORTRAN V for use on the 
UNIVAC 1108 under the EXEC II operating system. 

Documentation of the SEPS program is divided into two separate volumes. 

VOLUME I - Program Manual contained herein 

VOLUME II - User's Manual contains the information necessary for a 
user to adequately understand and use the SEPS Computer 
Program 
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1.0 INTRODUCTION 



1.0 INTRODUCTION 


This document contains information pertaining to the Program Manual, 
Programmer Guide, and Program Utilization of the Shuttle Electrical' 

Power System (SEPS) computer program. The main objective of this manual 
is to provide the information necessary to interpret and use the routines 
comprising the SEPS program. 

The subroutine descriptions are divided into four categories; control, 
Phase I, Phase II, and analysis routines. The subroutine descriptions 
include the name, -purpose, method (if applicable), variable definitions 
and logic flow. 

The SEPS User's Manual provides the information necessary for a 
user to adequately understand and use the SEPS computer program. 
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2.0 PROGRAM 

DESCRIPTION 



2.0 PROGRAM DESCRIPTION 


2.1 PURPOSE 

The TRW Shuttle EPS Analysis Computer Program (SEPS) was developed 
for use as a premission evaluation tool. The purposes of the program 
are to (a) predict EPS performance and EPS consumables usage when the 
system is operated and perturbed in accordance with premission flight 
plans, and (b) perform parametric and special case studies on the Shuttle 
EPS. 

2.2 SEPS PHASE I 
2.2.1 Discussion 

The SEPS Computer Program Phase I, using a mission event timeline, 
develops an electrical load profile and provides subsystem and mission analyses 
of the power and energy demands for Shuttle missions. The analysis is 
based on a 28 VDC load bus voltage. The. program utilizes a data base 
which describes all Shuttle electrical power consuming equipment in terms 
of power requirements and relating all the equipment to subsystems. This 
data combined with a desired mission event timeline provides the basis for 
the output interface tape consisting of event time point data and a listing 
of the activated components. The interface tape is utilized as the input 
driver for the Phase I COMUSE analyses and Phase II. The Phase I COMUSE 
analyses are discussed in more detail in the User's Manual, Section 2.1 
(Formatted Printout Description). 

The Phase I output and Phase I COMUSE output analyses have been adjusted 
to include an assumed 4% average line loss and an inverter efficiency of 80%. 
These data have been hard coded into the program and would require a change in 
Subroutine JVSEPS to revise these assumptions. The respective words are RESLOS 
for the line loss factor and PFEFF for the inverter efficiency. The line loss 
factor and inverter efficiency are not used in the Phase II analyses. 


p$B0)TNG * AGE BLANK NOP BBS® 
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2.3 PHASE II 
2.3.1 Discussion 


The SEPS Phase II program integrates the various math models that 
define the operating characteristics of the power sources, distribution 
and equipment of the Shuttle electrical power system. This provides the 
capability of simulating the total electrical power system with which 
system design and design/mission requirements compatibility can be analyzed 
and parametric studies performed. Another capability of SEPS Phase II is 
the fuel cell cryogenic requirements which result from the mission profile 
electrical power demand and operating procedures. 

The SEPS Phase II program utilizes the Phase I interface tape to 
provide the electrical load profile and active components. The user also 
has the capability to change or modify the configuration or input data 
through the use of an input card alter deck. 

2.3.2 Math Models 

In order to accomplish the SEPS Phase II capabilities several math 
models were required. A description of the math models, their intended 
use, and primary subroutine follows. 

EPS Distribution Circuit Math Model - This model is a representation 
of the Shuttle dc electrical power distribution and control system. 

Using a node analysis technique and the load profile contained on the 
Phase I interface tape, the distribution circuit model will determine 
the system bus voltages and currents and the load voltages and currents. 

The primary subroutine is DCSOLV. 

Fuel Cell Math Model - This model is a representation of the Shuttle 7 
KW average, 12 KW peak fuel cells. The fuel cell model is called by the 
distribution circuit model to provide source voltage as a function of 
source current. The fuel cell model also provides the cryogenic use 
rates to the cryogenics model. The primary subroutines are FUELIV and 
FUCLTM. 
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Inverter System Math Model - The inverter model is a simplified 
representation of the Shuttle 9 inverter/3 phase ac system. This model 
will calculate the inverter no load and load losses and reflect the ac 
inverter loads to the dc distribution system. The primary subroutine is 
ACINVT. 

Constraints Model - The constraints model will provide for automatic 
program checking and flagging of distribution, power source, and reactant 
storage system constraint violations. As actual performance, test and 
limit data becomes available, the constraints model can be updated. The 
primary subroutine is REDLIN. 

2.3.3 Phase I/Phase II Interface 

2.3.3. 1 Discussion 

The SEPS Phase I/Phase II interface has been designed to allow 
maximum utilization of analysis subroutines, output subroutines, and 
user interface through program option and control cards. Section 2. 3. 3. 2 
flow charts the Phase I/Phase II program. The flow chart shows a com- 
monality in program output and the use of the same analysis subroutines for 
both Phase I and Phase II. The user through program option and control 
cards can either allow or suppress virtually any portion of the combined 
Phase I/Phase II program outputs. The output interface tape of Phase I 
is the input driver for Phase II. Modifications to the interface tape 
data can be inserted by the user through an input card deck. 
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2. 3. 3. 2 Functional Flow Diagram 
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2. 3. 3. 2 Functional Flow Piaqram (Continued) 
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2.4 EXTERNAL PROGRAM INTERFACES 


To facilitate the SEPS program in meeting its defined requirements, 
the following interfacing programs have been developed: 

1. CIFLIS 

This program is used to list a given number of files from a 
card image tape. 

2. CMPDAT 

This program is used to give a component and/or mission phase 
comparison of two interface tapes 

3. CREDJR 

This program is used to make changes to the component definition 
card image tape file. 

4. FILDUP 

This program is used to duplicate card image files. 

5. JVMMPS 

This program is used to concentrate several timeline files 
together and to time order the resulting file. 

6. NEWHLP 

This program is used to construct a tape able to be plotted 
showing user designated component's time history of operation. 

7. STLPLT 

Generates CALCOMP plots from the unformatted interface output 
tapes. 

8. WLCCIT 

This program is used to create or alter a card image tape file, 
as a by-product the file is listed and each entry is numbered. 
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3.0 SUBROUTINE 
DESCRIPTION 



3.0 SUBROUTINE DESCRIPTION 


3.1 CONTROL SUBROUTINE 

3.1.1 Subroutine: JVSEPS 


PURPOSE: This. routine controls the SEPS program execution. 

METHOD: This routine controls the following functions: 

1. Determines total area of random access available 

2. Reads the option-units cards 

3. Reads the abort time 

4. If applicable, calls for the mission phase definition 
cards to be read 

5. If applicable, calls for the compacted component 
dictionary to be read for use in a Phase II only run 
utilizing an interface tape 

6. Controls the execution of 

a. Phase I 

b. Phase I COMUSE 

c. Phase II 

d. Phase II COMUSE 

.VARIABLES: The variables used in this subroutine are listed in the 

common blocks of the functional flowchart. Figure 3.1.1. 
See Appendix for definition of all variables. 


9 




CONT. ON PG 2 JVSEPS 

eg I. . qp— _z 

FIGURE 3 . 1 . 1 . FUNCTIONAL FLOWCHART OF SUBROUTINE: JVSEPS 


ORIGINAL PAGE IS 
OF POOR QUALET*! 
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C0NT. ON PG 3 JVSEPS 

PG ? OF 7 

.FIGURE 3.1.1 FUNCTIONAL FLOWCHART OF SUBROUTINE: JVSEPS (CONTINUED) 
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FORMAT C 1H1/40U0X. "IFLA6". I2.3X. 12. IX. 3A6/1/1 OX. "TDELT 1 ".FI 0.4/ 
10X. " TDELT A". 3X, FI 0.4/1 OX. "NPRT ",3X. t 10/1 OX. " JPR I NT” . 3X. 1 1 0/ 
1H1/////20U0X. I5.2X. A1 ,2X. "UNIT ", 15. 5X. “F ILE ‘. 15. 5X. 2A6/)/ ) 


IFCIFILEC 7D.6T.0_^ 



JFWLG25+IR^ 

I CALL MHANDL I 


CONT. ON PG 4 


JVSEPS 
PG 3 OF 


FIGURE 3.1.1. FUNCTIONAL FLOWCHART OF SUBROUTINE: JVSEPS (CONTINUED) 


ORIGINAL PAGE IS 
OP POOR QUALITY 












FIGURE 3 . 1 . 1 . FUNCTIONAL FLOWCHART OF SUBROUTINE: JVSEPS (CONTINUED) 
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CONT. ON PG G JVSEPS 

PG 5 OP 7 


FIGURE 3.1.1. FUNCTIONAL FLOWCHART OF SUBROUTINE: JVSEPS (CONTINUED) 
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CONT. ON PG 7 JVSEPS 

PG 6 OF 7 


FIGURE 3.1.1. FUNCTIONAL FLOWCHART OF SUBROUTINE: JVSEPS (CONTINUED) 
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FIGURE 3.1.1. FUNCTIONAL FLOWCHART OF SUBROUTINE: OVSEPS (CONTINUED) 
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3.2 PHASE I SUBROUTINES 
3.2.1 Subroutine: PHASE! 


PURPOSE: 

METHOD: 


VARIABLES: 


To control the creation of a load profile 

This routine controls the following functions 

1. Creates a component dictionary 

2. Creates a procedure dictionary 

3. Creates an activity dictionary 

4. Reads a timeline consisting of activities, procedures, 
components, switches, and cyclic elements and 'converts 
it to a component event timeline 

5. Compacts out the unused components in the component 
dictionary 

6. Analyzes the component event timeline 

The variables used in this subroutine are listed in the 

common blocks of the functional flowchart, Figure 3.2.1. 

See Appendix for definition of all variables. 
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SUBROUTINE PHASE 1 ( IHQR * NWL .IDA) 

INCLUDE STRAGA 

INCLUDE STRAG1 

INTEGER QQASG 

IU= 1UN I T( 1) 

REWIND IU 

CALL QQFILCIUNITC D’.IFILEC 1 3-1 . IS 


I 




CGNT, ON PG 2 

FIGURE 3.2,1 FUNCTIONAL FLOWCHART C 
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PHASE 1 

PR 1 OF 7 


SUBROUTINE PHASE I 









FIGURE 3 



CONT. ON PG 3 PHASE I 

P6 2 OF 7_ 

2.K. FUNCTIONAL FLOWCHART OF SUBROUTINE PHASE I (CONTINUED) 


ORIGINAL! PAGE IS 
OF POOR QUALM 
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FIGURE 



CONT. ON PG 4 


PHASE 1 

PG 3 OF 1 _ 


L2.1. FUNCTIONAL FLOWCHART OF SUBROUTINE PHASE I (CONTINUED) 
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CONT. ON PG 6 

FIGURE 3.2.1. FUNCTIONAL FLOWCHART OF SUBROUT 




of POOR QUALITY 



PHASE I (CONTINUED) 





C0NT. ON PG 7 


PHASE 1 


FIGURE 3.2.1. FUNCTIONAL FLOWCHART OF SUBROUTINE PHASE I (CONTINUED) 
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PHASE 1 


FIGURE 3.2.1 . FUNCTIONAL FLOWCHART OF SUBROUTINE PHASE I (CONTINUED) 


ORIGINAL PAGE IS 
OF POOR QUALM 
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3.2.2 Subroutine: ACYCLE 


PURPOSE: 

METHOD: 

VARIABLES: 

NOTE: 


This routine converts an Activity into Procedures and 
Components . 

This routine interrogates the Activity dictionary and 
calls the appropriate subroutines to correctly handle 
Procedures and Components. If the Activity cannot be 
found, the following diagnostic is generated. 

REQUESTED ACTIVITY NNNNNNNNNN IS NOT IN THE DICTIONARY 

The variables used in this subroutine are listed in the 
common blocks of the functional flowchart. Figure 3,2.2. 
See Appendix for definition of all variables. 

Subroutine ACYCLE is essentially identical to Subroutine 
AHANDL. The requirement for these subroutines is dictated 
by the program logic. 
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SUBROUTINE ACYCLECN. TT J 
INCLUDE STRAGA 
DIMENSION IA0ATC9) 

COMMON /UNITS/ IU5.IUG.IU7.IU8.IU9.IU10.IU11 

COMMON /ONE/ TYPA.TYPy.NUMA. MO. STRIA. STOPA. UFA. PERA.PONA 

equivalence itypa. iadatc m 

DATA NWAD /9/ 

IFQUND=0 



CONT. ON PG 2 ACYCLE 

P S 1 OF ; 

FIGURE 3.2,2. FUNCTIONAL FLOWCHART OF SUBROUTINE ACYCLE 
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ACYCLE 

P 6 2 FINAL 


FIGURE 3.2.2. FUNCTIONAL FLOWCHART OF SUBROUTINE ACYCLE (CONTINUED) 
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3.2.3 Subroutine: AHANDL 


PURPOSE: This routine converts an Activity into Procedures, 

Components and Switches. 

METHOD: This routine interrogates the Activity dictionary and 

calls the appropriate subroutines to correctly handle 
Procedures, Components, and Switches. If the Activity 
cannot be found, the following diagnostic is generated. 

REQUESTED ACTIVITY NNNNNNNNNN IS NOT IN THE DICTIONARY 

VARIABLES: The variables used in this subroutine are listed in the 

common blocks of the functional flowchart. Figure 3.2.3. 
See Appendix for definition of all variables. 


PRECEDING PAGE BLANK NOT FTT.M^n 
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SUBROUTINE AHANOLCN.TT] 

INCLUDE STRAGA 
DIMENSION IADATO) 

COMMON /UNITS/ IUS.IUE.IU7.IU8.IU9.IU10.IU11 
COMMON /ONE/ TYPA.TYPr.NUMA.MO.STRTA.STOPA.UFA.PERA.PONA 
EQUIVALENCE C TYPA. IADATC 1)3 
DATA NWAD /9/ 

IFOUNO=0 



CONT. ON PG 2 


AHANDL 


FIGURE 3.2.3. FUNCTIONAL FLOWCHART OF. SUBROUTINE AHANDL 


ORIGINAL PAGE IS 
•OF POOR QUALITY 
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100 TO 15" 



CONT. ON PG 3 AHANDL 

P6 2 . S£ 4 

FIGURE 3.2.3. FUNCTIONAL FLOWCHART OF SUBROUTINE AHANDL (CONTINUED) . 


DBIGETAE PAGE IS 

<W POOR QUALITY 
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CALL CHANOLCN. 0.0. .0.. 0. . NUMA. MD. STRT A. STOPA. UFA. 
TT) 



FIGURE 3.2.3. 


CONT. ON P6 4 


AHANDL 


FUNCTIONAL FLOWCHART OF SUBROUTINE AHANDL (CONTINUED) 


prightae page 2B 

op POOR QUALITY 
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AHANDL 


FIGURE 3.2.3. FUNCTIONAL FLOWCHART OF SUBROUTINE AHANDL (CONTINUED) 
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3.2.4 Subroutine: AREAD 


PURPOSE: This routine reads the Activity/Procedure definition 

cards and creates the Activity/ Procedure dictionary. 

METHOD: This routine is called twice; once to create the Activity 

dictionary and once to create the Procedure dictionary. 

Each card as it is read in is tested to determine if it 
is a File Title card, a Card Type 1 (an Activity/ Procedure 
title card), or a Card Type 2 (an Activity/Procedure 
definition card). As each Card Type 2 is read, it is 
related to its card type 1 and stored in a Sequentially/ 
indexed file and its location in the file is stored in 
the dictionary. 

The following diagnostics are provided: 

INCORRECT TYPE AANNNNNNN NOT 

NO MORE DRUM SPACE AVAILABLE - THE LAST CARD PROCESSED WAS AA 


ID NNNNNNN NUMBER NNNNNNNNNN RUN ABORTED 
READ ERROR ON CARD AANNNNNNN 
VARIABLES: Calling Arguments: 


IDIOT - where the dictionary is to be stored 
ICD - number of entries in IDICT 
NWL - words of drum remaining 

IDA - drum address 

IU - unit definitions' are to be read form 

ITY - entry type 1 - procedure 2 - activity 

IDEM - IDICT dimension 


The remaining variables listed in this subroutine are 
listed in the common blocks of the functional flowchart. 
Figure 3.2.4. See Appendix for definition of all variables. 




£S> 
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READUU. 1000. END= 100. ERR=30 D IT* IDY. ID.MD. I STRT . ISTQP.UFS. 


PERS.PGNS. I SWT 


VR1TE( IU6.2000) IT. IDY. ID.MD. ISTRT . 1STOP. UFS. PERS. PONS. ISWT 



V 


v 

C0NT. ON PG 2 AREAD 

1 PG 1 OF a 


FIGURE 3,2.4. FUNCTIONAL FLOWCHART OF SUBROUTINE AREAD 
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CONT. ON P 6 3 


ARE AD 


FIGURE 3.2.4. FUNCTIONAL FLOWCHART OF SUBROUTINE AREAD (CONTINUED) 


ORIGIN AD PAGB IS 
OF POOR QUALITY 
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FIGURE 3.2.4. FUNCTIONAL FLOWCHART OF SUBROUTINE AREAD (CONTINUED) 
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CQNT. ON PG 5 ARE AO 

PG 4 OP 8 


FIGURE 3.2.4. FUNCTIONAL FLOWCHART OF SUBROUTINE AREAD (CONTINUED) 
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1020 j. 

FORMAT (F5. 1 D 

^ 

UF=UF/1 00. 

DEC0DEC6. 1025.PERS) PER 
DECODECS. 1025. PONS) PON 


1025 & 
FORMAT CF6. 0) 


IDC=!OC+ 1 
ID1CTC IOC. 1)=IAID 
IOICTC IDC. 2)=I0A 


CALL RWRITEC IOA. IDAT.NWAD. ISTAT) 



CONT. ON PG S ARE AD 

PS_ 5 , OF S_ 

FIGURE 3.2.4. FUNCTIONAL FLOWCHART OF SUBROUTINE AREAD (CONTINUED) 




Of 1 
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FIGURE 3.2.4. FUNCTIONAL FLOWCHART OF SUBROUTINE AREAD (CONTINUED) 
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CONT. ON PG 8 


AREAO 


FIGURE 3.2.4. FUNCTIONAL FLOWCHART OF SUBROUTINE AREAD (CONTINUED) 
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FIGURE 3.2.4. FUNCTIONAL FLOWCHART OF SUBROUTINE AREAD (CONTINUED) 
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3.2.5 Subroutine: CCYCLE 


PURPOSE: This routine creates the component portion of the event 

timeline. 

METHOD: For each component the following are determined: 

1. Location in the component dictionary 

2. Shows the component as "active" and to be included in 
the compacted dictionary 

3. Determines the event on and off time 

4. Stores the event in the timeline array 
- 5. Writes the event on drum 

VARIABLES: The variables used in this subroutine are listed in the 

common blocks of the functional flowchart. Figure 3.2.5. 
See Appendix for definition of all variables. 


NOTE: Subroutine CCYCLE is essentially identical to Subroutine 

CHANDL. The requirement for these subroutines is dictated 
by the program logic. 


preceding 


page BLANK HOESJa*® 


- 
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SUBROUTINE CCYCLECN. NA. SA.EA.UA. NP. HP. SP. EP. UP. TT ) 
INCLUDE STREGA 

COMMON /UNITS/ IU5. IUG. IU7. IU8. IU9. IUtO, IU1 I 
COMMON /TL INF/ I OUT. IOUTM, IDRM. I END. IF IL. IDA.NWL. ITL 
DIMENSION ITL(5000.2) 

DIMENSION TTLC5000.2) 

EQUIVALENCE! ITLC 1. 1 3.TTLCI. 1 ) ) 

IA=1 



CONT. ON P6 2 


CCYCLE 


FIGURE 3.2.5. FUNCTIONAL FLOWCHART OF SUBROUTINE CCYCLE 
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SO to 50 >Hlz 

< 



V/RlTEt IUG. 9000D NP.NA.N.TT 

1 


9000 i 

1 


| FORMAT (5X* ,l COMPONENT “•110*“ IS NOT IN THE D I CT IONARY M •2115.F10.4) 


I 60 T0~5Q~> ^^ 

20 i 


CONTINUE 


CALL TIHSETCN.TT.SA.EA.SP.EP.S.Ei 


H 


ICUT=IGUT+1 
TTLC I OUT . I )=S 
I TLC I OUT . 2)= IDA 
NS=1 


±j 


IFCN . NE. 0) 


J 

T 

z 

| N$=NS*N / 1 ABSCN ) 

T 

7 

*5 


CONT. ON P6 3 


CCYCLE 

PG 2 OP 5 


FIGURE 3.2.5. FUNCTIONAL FLOWCHART OF SUBROUTINE CCYCLE (CONTINUED) 


' page 18 
Quality 
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CONT. ON PG 4 


CCYCLE 


FIGURE 3.2.5. FUNCTIONAL FLOWCHART OF SUBROUTINE CCYCLE (CONTINUED) 
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CALL DRMRITCNWL. 




ORIGINAL PAGBJ3' 
OF POOR QUAUT2 



IAL FLOWCHART OF SUBROUTINE CGYCLE (CONTINUED) 


i 
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CCYCLE 


FIGURE 3.2.5. FUNCTIONAL FLOWCHART OF SUBROUTINE CCYCLE (CONTINUED) 
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3.2.6 Subroutine: CHANDL 


PURPOSE: 

METHOD: 


VARIABLES: 


This routine creates the component portion of the event 
timeline. 

For each component the following are determined: 

1. Location in the component dictionary 

2. Shows the component as "active" and to be included in 
the compacted dictionary 

3. Determines the event on and off time 

4. Stores the event in the timeline array 

5. Writes the event on drum 

The variables used in this subroutine are listed in the 
common blocks of the functional flowchart. Figure 3.2.6. 
See Appendix for definition of all variables. 
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SUBROUTINE CNANDLCN.NA.SA.EA.UA.NP.MP.SP.EP.UP. TTJ 
INCLUDE STRAGA 

COMMON /UNITS/ IU5. IUS. IU7. IU8. IU9. IU1 0. IU1 I 
COMMON /TL INF/ I OUT . I OuTM. 1 DRM. I END. IF I L . I DA. NWL . I TL 
COMMON /FCYCL/ KF 

COMMON /ALTERN/ ICDCBC 750) . I CDCAC500. 3] . I ACT 
DIMENSION ITU50Q0.2) 

DIMENSION TTU5000. 2) 


K 



FIGURE 3.2.6. FUNCTIONAL FLOWCHART OF SUBROUTINE CHANDL 
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FIGURE 3.2.6. FUNCTIONAL FLOWCHART OF SUBROUTINE CHANDL (CONTINUED) 
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FIGURE 3.2.6. FUNCTIONAL FLOWCHART OF SUBROUTINE CHANDL (CONTINUED) 



54 









FIGURE 3.2.6. FUNCTIONAL FLOWCHART OF SUBROUTINE CHANDL (CONTINUED) 
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3.2.7 Subroutine : COMPCT 


PURPOSE: 

METHOD: 

VARIABLES: 


To eliminate unused components from the component 
di cti onary 

Only components found to be “active" are stored in 
the compacted dictionary. If requested, the compacted 
dictionary is written to an auxiliary unit for use in 
a later execution. 

The variables used in this subroutine are listed in the 
common blocks of the functional flowchart. Figure 
3.2.7. See Appendix for definition of all variables. 
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SUBROUTINE COMPCTflUJ 
INCLUOE STRAGA 
INTEGER OOASG 
DIMENSION IDATC50) 

COMMON /ALTERN/ ICDCBC75CO . ICDCAC500. 3) . I ACT 
COMMON /BSLOCA/ IOUMBC750). IDUMSC750) 

COMMON /BSLOC / IBUSCC500) . ISYSCC500) 



FIGURE 3.2.7. FUNCTIONAL FLOWCHART OF SUBROUTINE COMPCT 
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FIGURE 3. 



CQNT. ON PG 3 CGMPCT 

PQ 2 0 E 2 

.7. FUNCTIONAL FLOWCHART OF SUBROUTINE COMPCT (CONTINUED) 
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catiPCT 

PG 2_ F 1NA1. 


FIGURE 3.2.7. FUNCTIONAL FLOWCHART OF SUBROUTINE COMPCT (CONTINUED) 


ORIGINAL PAGE IS 
OF POOR' QUALITY? 
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3.2.8. Subroutine: CREAD 


PURPOSE: 

METHOD: 

VARIABLES: 


This routine creates the component dictionary. 

The component definition is read. The component modes 
are gathered together and written randomly on drum. 

The component dictionary consists of component ID number 
and drum location. 

The variables used in this subroutine are listed in the 
common blocks of the functional flowchart. Figure 3.2.8. 
See Appendix for definition of all variables. 
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FGRMAT(5X. A 1 . 11. IX. 16. I2.2X.5A6. IX.F5. 0.4F7. 0. 13) 


CGNT. ON PG 2 


CREAD 


FIGURE 3.2.8. FUNCTIONAL FLOWCHART OF SUBROUTINE CREAD 
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IFCIFRST.LT.l) 


F 


T 



FIGURE 3.2.8. FUNCTIONAL FLOWCHART OF SUBROUTINE CREAD (CONTINUED) 
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I IFRST=IFRST+1 

ILNCP= ILN* 1 000+ ICP 
IBUSCC IDC+n = lLN 
ISYSCC IDC+1)=ID/10000 
NAMS(U=NAMCn 
NAMS C 2 ) =NAM ( 2 ) 
NAMSC3)=NAMC3) 
NAMSt4)=NAM£ 4 ) 
NAMSC5J=NAM(5) 


NAMSC6)=6H 

PFS=PF 


V 

L BO .J Q . 
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CQNT. ON PG 4 


CREAO 

PG 3 OF 6 


FIGURE 3.2.8. FUNCTIONAL FLOWCHART OF SUBROUTINE CREAD (CONTINUED) 


QRIGINAE PAGE m 
OF POOR QUALITY! 
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CONI. ON PG 5 


CREAO 


FIGURE 3.2.8. FUNCTIONAL FLOWCHART OF SUBROUTINE CREAD (CONTINUED) 


65 












, 5070 ± 

FGRMATC5X. "READ ERROR ON CARD ". A2, 17. 5X. "COMPONENT" ) 

_ T _ 

CONT. ON PG G CREAD 

P6 .. 5 . . 0 F 6 

FIGURE 3.2.8. FUNCTIONAL FLOWCHART OF SUBROUTINE CREAD (CONTINUED) 
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END 


'ORIGINAL PAGE IS 
OF POOR QUALM 


CREAO 

PS-.6 PINAL 


FIGURE 3.2.8. FUNCTIONAL FLOWCHART OF SUBROUTINE CREAD (CONTINUED) 







INTENTIONALLY 

LEFT 


BLANK 



3.2.9 Subroutine: CTAPE 


PURPOSE: 

METHOD: 


VARIABLES: 


To write the Phase I interface tape. 

The changes of component loads and switch positions 
from the previous time are determined. Only the changes 
are written on the interface tape. The first record 
contains zero values for all switch positions. The 
first record is used for initialization. 

The variables used in this subroutine are listed in the 
common blocks of the functional flowchart. Figure 2.3.9. 
See Appendix for definition of all variables. 
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CONT. ON PG 2 CTAPE 


FIGURE 3.2.9. FUNCTIONAL FLOWCHART OF SUBROUTINE CTAPE 


IrbdedWg WlGB blank NOS 0WSED! 


70 










CONT. ON PG 3 


CTAPE 


FIGURE 3.2.9. FUNCTIONAL FLOWCHART OF SUBROUTINE CTAPE (CONTINUED) 


*SS8S? 
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CGNT. ON PG 4 


CTAPE 


FIGURE 3.2.9. FUNCTIONAL FLOWCHART OF SUBROUTINE CTAPE (CONTINUED) 
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FIGURE 3.2.9. FUNCTIONAL FLOWCHART OF SUBROUTINE CTAPE (CONTINUED) 
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WRITECIU) CTIME.PF.NW. CNTAPEC ID . C TAPELU . J D . J= 1 . 3) . l=l.NW) 
NWS, CNTAPSC I). ITAPStn. 1 = 1. NWS) 



D 
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3.2.10 Subroutine: CYCLE 


PURPOSE: 

METHOD: 

VARIABLES: 


To control all cyclic elements. 

When a cyclic element is encountered it is stored in the 
cyclic definition table. Periodically this routine is 
called to update the cyclic's condition by calling either 
AHANDL, PHANDL, or CHANDL depending upon the type of 
cyclic element. 

The variables used in this subroutine are listed in the 
common blocks of the functional flowchart. Figure 3.2.10. 
See Appendix for definition of all variables. 
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SUBROUTINE CYCLE! T.N. MO, PER. PON. TT.KIND) 
DIMENSION TSCIOO) 

DIMENSION NSC I 00) 

DIMENSION MSC100) 

DIMENSION PERSUOO) 

DIMENSION PONS C 1 00) 

DIMENSION TTS(IOO) 

DIMENSION NSTORCIOO.S) 


COMMON /TL INF / IOUT. IOUTM 
COMMON /CYCLIC/ NS. MS. TS. PERS. PONS. TTS. I 
COMMON /UNITS/ IU5. I US. IU7. IU8. IU9. IU1 0. lull 
COMMON /FCYCL/ KF 
DATA I /0/ 

DATA NS /l 00*9999999/ 

EQUIVALENCE CNSC 1) . NSTORC 1 . 1) ) 



FIGURE 3.2.10. FUNCTIONAL FLOWCHART OF SUBROUTINE CYCLE 
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PERSC I )=PER 
PONS(n=PON 
TTSCI)=TT 

" "l '• 

Treturn I 



CONT. ON PG 3 CYCLE 

P6 2 OF & 


FIGURE 3.2.10. FUNCTIONAL FLOWCHART OF SUBROUTINE CYCLE (CONTINUED) 









FIGURE 3.2.10. FUNCTIONAL FLOWCHART OF SUBROUTINE CYCLE (CONTINUED) 
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C0NT. ON P G 5 


CYCLE 

PG 4 OF £ 


FIGURE 3.2.10. FUNCTIONAL FLOWCHART OF SUBROUTINE CYCLE (CONTINUED) 
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CONT. ON PG G 


CYCLE 


FIGURE 3.2.10. FUNCTIONAL FLOWCHART OF SUBROUTINE CYCLE (CONTINUED) 


8 








CALL CHANDLCL. 0. 0. . 0. . 0. .NSCK). MSCKJ . 0. . 0. . 0 
TTSCK) ) 



FIGURE 3.2.10. FUNCTIONAL FLOWCHART OF SUBROUTINE CYCLE (CONTINUED) 


0RIGINAI1 PA.G3HS 
OF POOR QUAIOT 


81 











3.2.11. Subroutine: MHANDL ' 

PURPOSE: Create the Mission Phase Definition Dictionary 

METHOD: Read the Mission Phase Definition cards and store 

them in an array for later use. 

VARIABLES: The variables used in this subroutine are listed in the 

common blocks of the functional flowchart, Figure 3.2.11. 
See Appendix for definition of all variables. 


FBSC0MNG PAGE BLANK N03! 
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SUBROUTINE MHANOL 
INCLUDE STRA6A 

COMMON /UNITS/ IU5. IU6. IU7 . IU8 , I U3. 1U1 0. IU 1 t 
COMMON /MTRANT/ T IMEMC50) 

DIMENSION T l TL ( 6 ) 

DATA TiTLM /3 00*6H / 

DATA TM /50*9999339999./ 

MM-0 



L. 

AO \ J 2 — ^ — 

2 

| CONTINUE | 


READC IU5. 1 000. EN0=4 0< ERR=30) TYPE, OT. TIME. T I TL 

1 000 ^ 




IFCTYPE,NE,2HMM) 



2 

1 GO TO 20 ^ 



l 

2 

| MM=MM+ l 


Xffl 


TM(MM)=TIME 
TIMEMtMM)=OT 
TITLHU.MMD=TlTL(n 
T ITLMC 2. MM)=T [ TL( 2) 
TJTLMC3. MM) =T I TLC 3} 
T ITLMC 4. MM) =T ITL ( 4 ) 
T I TLMt 5.MM) = TITL(5) 
TITLMCS.MM)=TITLCS) 


T 


CONT. 


ON PG 


MHANOL 
PQ_J QE_ 


FIGURE 3.2.11. FUNCTIONAL FLOWCHART OF SUBROUTINE MHANDL 
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CGNT. ON PG 3 MHANDL 

PC? .2.-QE 3 

FIGURE 3.2.11. FUNCTIONAL FLOWCHART OF SUBROUTINE MHANDL (CONTINUED) 
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MHANDL 

PG. 3 F.IJjAL 


FIGURE 3.2.11. FUNCTIONAL FLOWCHART OF SUBROUTINE MHANDL (CONTINUED) 
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3.2.12 Subroutine: PCYCLE 


PURPOSE: 

METHOD: 

VARIABLES: 

NOTE: 


This routine converts Procedures into Components. 

This routine interrogates the Procedure dictionary and 
calls the appropriate subroutine to correctly handle 
Components. If the Procedure cannot be found, the 
following diagnostic is generated. 

REQUESTED PROCEDURE NNNNNNNNNN IS NOT IN THE DICTIONARY 

The variables used in this subroutine are listed in the 
conmon blocks of the functional flowchart. Figure 3.2.12. 
See Appendix for definition of all variables. 

Subroutine PCYCLE is essentially identical to Subroutine 
PHANDL. The requirement for these subroutines is dictated 
by the program logic. 
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SUBROUTINE PCYCLECN.NA.S.E.U.TT) 

INCLUDE STRAGA 
DIMENSION IPDATC9) 

COMMON /UNITS/ IU5. IU6. IU7. IU8. IU9. IUI0. IUH 

COMMON /TWO/ TYPP.TYPY.NUMP.MP.STRTP.STOPP.UFP.PERP.PONP 

EQUIVALENCE C TYPP, IPDATC 1 ) ) 

DATA NWPD /9/ 

IFOUND=0 



FIGURE 3.2.12. FUNCTIONAL FLOWCHART OF SUBROUTINE PCYCLE 
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j END | 


PCYCLE 

PG ? FINAL 


FIGURE 3.2.12. FUNCTIONAL FLOWCHART OF SUBROUTINE PCYCLE (CONTINUED) 
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3.2.13 Subroutine: PHANDL 


PURPOSE: This routine converts a Procedure into Components and 

switches. 

METHOD: This routine interrogates the Procedure dictionary and 

calls the appropriate subroutines to correctly handle 
Components and Switches. If the Procedure cannot be 
found, the following diagnostic is generated. 

REQUESTED PROCEDURE NNNNNNNNNN IS NOT IN THE DICTIONARY 

VARIABLES: The variables used in this subroutine are listed in the 

common blocks of the functional flowchart. Figure 3.2.13. 
See Appendix for definition of all variables. 




G Page 


Blank 
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SUBROUTINE PHANDLCN.NA.S.E.U.TT) 

INCLUOE StRAGA 
DIMENSION IPDATC9) 

COMMON /UNITS/ IU5.IUS.IU7.IU8.IU9.IU10.IU11 
COMMON /TWO/ TYPP.TYPY.NUMP.MP.STRTP.STOPP.UFP.PERP.PONP 
EQUIVALENCE CTYPP. IPDATt m 
DATA NWPD /9/ 

1FOUND=0 



CONT. ON PG 2 


PHANDL . 

PQ_J Q£ S_ 


FIGURE 3.2.13. FUNCTIONAL FLOWCHART OF SUBROUTINE PHANDL 

\ 


ORIGINAL 
OF POOR QTJAIiM 
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FIGURE 3.2.13. FUNCTIONAL FLOWCHART OF SUBROUTINE PHANDL (CONTINUED) 
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FIGURE 3.2.13. FUNCTIONAL FLOWCHART OF SUBROUTINE PHANDL (CONTINUED) 
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FIGURE 3.2.13. FUNCTIONAL FLOWCHART FOR SUBROUTINE PHANDL (CONTINUED) 


ORIGINAL M5GBH 
OF POOR QUALM. 
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PHANDL 


FIGURE 3.2.13. FUNCTIONAL FLOWCHART OF SUBROUTINE PHANDL (CONTINUED) 
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3.2.14 Subroutine: PLTNOW 


PURPOSE: Write the plot of total source power versus time, 

as determined in Phase I, on the printer. 

METHOD: Create headings necessary for the plot and write the 

plot on the printer. 

VARIABLES: The variables used in this subroutine are listed in the 

common blocks of the functional flowchart. Figure 3.2.14. 
See Appendix for definition of all variables. 


C -3, 
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subroutine PLTNOWCLPAG. 1U, TEND) 

DIMENSION LPAGC 111.51) 



CONT. ON PG 2 PLTNOW 


FIGURE 3,2.14. FUNCTIONAL FLOWCHART OF SUBROUTINE PLTNOW 


OBIGWAEPAGSE 

OF POOR. QUAIira 
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CONT. ON PG 3 


PLTNOW 


FIGURE 3.2.14. FUNCTIONAL FLOWCHART OF SUBROUTINE PLTNOW (CONTINUED) 
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C0NT. ON PG 4 PLTNOW 

PG 3 OF 5 

FIGURE 3.2,14. FUNCTIONAL FLOWCHART OF SUBROUTINE PLTNOW (CONTINUED) 
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.1020 j 
| FORMAT Cl HI) 


CfW T . ON PG 5 PLTNOW 

BS-4. OF 5 

FIGURE 3.2.14. FUNCTIONAL FLOWCHART OF SUBROUTINE PLTNOW (CONTINUED) 
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3.2.15 Subroutine: PRTPLT 


PURPOSE: 

METHOD: 

VARIABLES: 


Create a plot of total source power versus time, as 
determined in Phase I. 

Determine the scaling necessary for this point. Store 
the point in the plot array. If the time elapsed has 
exceeded 24 hours, call PLTNOW to print the plot. 

The variables used in this subroutine are listed in 

the common blocks of the functional flowchart. Figure 3.2.15. 

See Appendix for definition of all variables. 
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subroutine prtpltct.p. IU.K) 
oimension lpagec i u.sn 


OATA I FRST /Q/ 

DATA IYS /l/ 

DATA LPAGE /5G61*6H / 




CONT. ON PG 2 


FIGURE 3.2.15. FUNCTIONAL FLOWCHART 


PRTPLT 
PG 1 OF 


SUBROUTINE PRTPLT 









5 






FIGURE 3.2.15. FUNCTIONAL FLOWCHART OF SUBROUTINE PRTPLT (CONTINUED) 
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3.2.16 Subroutine: RCYCLE 


PURPOSE: 

METHOD: 

VARIABLES: 


To control all cyclic elements. 

Periodically this routine is called to update the 
cyclic's condition by calling either AHANDL, PHANDL, 
or CHANDL depending upon the type of cyclic element. 

The variables used in this subroutine are listed in the 
common blocks of the functional flowchart. Figure 3.2.16. 
See Appendix for definition of all variables 
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CONT. ON PC 2 


BCYCLE 


FIGURE 3.2.16. FUNCTIONAL FLOWCHART OF SUBROUTINE RCYCLE 
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FIGURE 3.2.16. FUNCTIONAL FLOWCHART OF SUBROUTINE RCYCLE (CONTINUED) 
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CQNT. ON PG 4 RCYCLE 

PG 3 OF 5 

FIGURE 3.2,16. FUNCTIONAL FLOWCHART OF SUBROUTINE RCYCLE (CONTINUED 






no 








CONT. ON PG 5 


RCYCLE 


FIGURE 3.2.16. FUNCTIONAL FLOWCHART OF SUBROUTINE RCYCLE (CONTINUED) 
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RCYCLE 

PG_5 FINAL 


FIGURE 3.2.16. FUNCTIONAL FLOWCHART OF SUBROUTINE RCYCLE (CONTINUED) 


ORIGINAL PAGE 
OF POOR QUALITY 
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3.2.17 Subroutine: RFIL 


PURPOSE: 

METHOD: 

VARIABLES: 


To handle the situation created when the event array 
overloads. 

If necessary, all cyclic elements status are brought up 
to the present time. The event array is ordered on 
time. The number of points to be written to the event 
timeline unit (0) is determined. The points are written 
to the tape and removed from the event array. 

The variables used in this subroutine are listed in the 
common blocks of the functional flowchart. Figure 3.2.17. 
See Appendix for definition of all variables. 
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[QUTM=5000 

< DO 10 1=1. NT > 

| CAUL RCYCLE C TSC i) . NSU) . T ImY” 

a 




CONTINUE 


0 5 


f 

CONTINUE 


} CALL ROROERCT. 5000. 5000. 2) 


=a 


CONT. ON PG 2 


RFIL 

P£_J QE 1 


FIGURE 3.2.17. FUNCTIONAL FLOWCHART OF SUBROUTINE RFIL 
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C0NT. ON PG 3 


RFIL 

PG 2 — OF 3L 


FIGURE 3.2.17. FUNCTIONAL FLOWCHART OF SUBROUTINE RFIL (CONTINUED) 
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RFIL 

ee_2 final 


FIGURE 3.2,17. FUNCTIONAL FLOWCHART OF SUBROUTINE RFIL (CONTINUED) 
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3.2.18 Subroutine: SHANDL 


PURPOSE: 

METHOD: 

VARIABLES: 


This routine creates the switch portion of the event 
timeline. 

For each switch the following are determined: 

1. Determines the event on and off time 

2. Stores the event in the timeline array 

3. Writes the event on drum 

The variables used in this subroutine are listed in the 
common blocks of the functional flowchart, Figure 3.2.18. 
See Appendix for definition of all variables. 
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FIGURE 3.2.18. FUNCTIONAL FLOWCHART OF SUBROUTINE SHANDL 
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CONT. ON PG 3 SHANOL 

PG 1 OF 5 

FIGURE 3.2.18. FUNCTIONAL FLOWCHART OF SUBROUTINE SHANDL (CONTINUED) 
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IQUT=IQUT+1 


TTLUOUT. 1)=E 
$ 

CONT. ON PG 4 SHANDL 

PG 3 OP 5 

FIGURE 3.2.18. FUNCTIONAL FLOWCHART OF SUBROUTINE SHANDL (CONTINUED) 


OBIGINAIi PAG® TS 

OF POOR QUAIiTTY! 
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FIGURE 3.?. 18. FUNCTIONAL FLOWCHART OF SUBROUTINE SHANDL (CONTINUED) 
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SHANOL 


FIGURE 3.2.18. FUNCTIONAL FLOWCHART OF SUBROUTINE SHANDL (CONTINUED) 


ORIGINAL! PAGETIgf 
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3.2.19 Subroutine: TPOUTJ 


PURPOSE: 

METHOD: 

VARIABLES: 


Analyze the event timeline and provide the interface 
tape, plots, and all formatted printouts from Phase I. 

The event timeline is read to determine the change of 
status for switches and components. If the change of 
status effects a component certain checks are made to 
insure the component is operating at its highest use 
factor. Changes in status happening within an input 
delta time of each other are grouped together and printed 
and plotted as if they occurred at the same time. 

The variables used in this subroutine are listed in the 
common blocks of the functional flowchart. Figure 3.2.19. 
See Appendix for definition of all variables. 


123 




CONT. ON PG 2 TPOUTJ 

Pfi 1 OF 22 

FIGURE 3.2.19. FUNCTIONAL FLOWCHART OF SUBROUTINE TPOUTJ 


ORIGINAL PAGE IS 
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FIGURE 3.2.19. 


CONT. ON PG 


TPOUTJ 


FUNCTIONAL FLOWCHART OF SUBROUTINE TPOUTJ (CONTINUED) 
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FIGURE 3.2.19. FUNCTIONAL FLOWCHART OF SUBROUTINE TPOUTJ (CONTINUED) 
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FIGURE 3.2. 



PWRC 1 J=PWRM( 1 .MD+ I ) 
PWRC2 }=PWRM( 2. M0 + I ) 
PWRC3)=PWRMC3.MD+ 1 ) 
UF=UST 



CGNT. GN PG 5 TPGUTJ 

P6 4 OF 32 

19. FUNCTIONAL FLOWCHART OF SUBROUTINE TPOUTJ (CONTINUED) 


127 






128 



SUBROUTINE TPOUTJ (CONTINUED) 












FIGURE 3.2.19. FUNCTIONAL FLOWCHART OF SUBROUTINE TPOUTJ (CONTINUED) 
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5 



CCNT. ON PG 8 


TPOUTJ 


FIGURE 3.2.19. FUNCTIONAL FLOWCHART OF SUBROUTINE TPOUTJ (CONTINUED) 
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FIGURE 3.2.19. FUNCTIONAL FLOWCHART OF SUBROUTINE TPOUTJ (CONTINUED) 
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FORMAT C5X. "***** COMPONENT ".3110." TURNED OFF INCORRECTLY " D 





CONT. ON PG 10 TPOUTJ 

P C , - 9- QF 22 

FIGURE 3.2.19. FUNCTIONAL FLOWCHART OF SUBROUTINE TPOUTJ (CONTINUED) 
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CONT. ON PG 11 TPOUTJ 

PG 10 OF 32 

FIGURE 3.2.19. FUNCTIONAL FLOWCHART OF SUBROUTINE TPOUTJ (CONTINUED) 
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} 




CONT. ON PG 12 TPOUTJ 

P6.,U QE 22 

FIGURE 3.2.19. FUNCTIONAL FLOWCHART OF SUBROUTINE TPOUTJ (CONTINUED) 
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CONT. ON PG 13 


TPOUTJ 


FIGURE 3.2.19. FUNCTIONAL FLOWCHART OF SUBROUTINE TPOUTJ (CONTINUED) 
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t> 15 


C0NT. ON PS 14 TPGUTJ 

PG 13 OF 22 

4 

FIGURE 3.2.19. FUNCTIONAL FLOWCHART OF SUBROUTINE TPOUTU (CONTINUED) 
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CGNT. ON PG 15 


TPOUTJ 


FIGURE 3.2.19. FUNCTIONAL FLOWCHART OF SUBROUTINE TPOUTJ (CONTINUED) 
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CGNT. ON PG 1G • TP0UTJ 

PG 1 5 Q£ 22 


FIGURE 3.2.19. FUNCTIONAL FLOWCHART OF SUBROUTINE TPOUTJ (CONTINUED) 


138 












C0NT. ON PG 17 TPOUTJ 

PS-16,. OF 22 


FIGURE 3.2.19. FUNCTIONAL FLOWCHART OF SUBROUTINE TPOUTJ (CONTINUED) 
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FIGURE 3.2.19. FUNCTIONAL FLOWCHART OF SUBROUTINE TPOUTJ (CONTINUED) 


ORIGINAL PAGE 16 
OF POOR QUALITY '' 
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im 


co 


1 


BUSCK3=BUS(K3+TLGADC J 3 


1 ? - 


5 60 


CONTINUE 


BUSC 2 3=BUS( 53+BUSt 83+BUSC 93+BUSC203 

BUSC 33=BUSC 63+BUSC 103+BUSC 1 1 )+BUS(2i 1 

BUSC 43=BUSC 73+BUSt 123+BUSC 133+BUSC223-t-BUSC233 

FC1=BUSC23 

FC2=8US(33 

FC3=BUSC43 

CRYQ=FC1 +FC2+FC3 

SLGAO=0. 


' < DO 57 0 "j = 1 .30 > 
^ 


SLOAD=SLOAO+TLOAD( J) 


5 70 


CONT I NUE 


~<Ldq sai) Jr . i .u2 > 
2 


JJ=JLGCACJ3 

PFBeJ3=TLGADCJJ3/SQRTCTLGAC3(JJ3*TLOAOCJJ3+ACLGAOCJ3*ACLGAOCJ33 
PFBC J3=S IGNCPFBC J 3 . ACLOADC J 3 3 



CONT. ON PG 20 TPOUTJ 

PG .12-QF 22 

FIGURE 3.2.19. FUNCTIONAL FLOWCHART OF SUBROUTINE TPOUTJ (CONTINUED) 
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FIGURE 



, 2 

CALL PRTPLTCOTIME.CRYO. IU6. 0) 



CALL CTAPECDT IME. TDELT1 . TR ILO. ICN. IU.PF8. I S 5 



WR I TE C IUP D DT I ME . SLQAO. CRY0 . FC 1 . FC2 . FC3. SUBSYS 




2Z 


3.2.19. FUNCTIONAL FLOWCHART OF SUBROUTINE TPOUTJ (CONTINUED) 
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FIGURE 3. 



[ WRITECIU6. 121 0) ( J, TLOADC J ) . J ~ 1 . 30) 



dJLmI>E 


< DO 605 J= 1 » 2 0 . 4 > - - - t> 22 
_s_ 


Jl=ILOCC JD 
J2=IL0C£J+1D 


TPOUTJ 

P(? 21 . Q F 32. 


CONT. ON PG 22 

.19. FUNCTIONAL FLOWCHART OF SUBROUTINE TPOUTU (CONTINUED) 
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C0NT. ON PG 23 TROUT J 

PS 22 J2E 22 

FIGURE 3.2.19. FUNCTIONAL FLOWCHART OF SUBROUTINE TPOUTJ (CONTINUED) 
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CONT. ON PG 24 


TPOUTJ 


FIGURE 3.2.19. FUNCTIONAL FLOWCHART OF SUBROUTINE TPOUTJ (CONTINUED) 
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CONT. ON PG 26 TPOUTJ 

P£_25 PF 22 

FIGURE 3.2.19. FUNCTIONAL FLOWCHART OF SUBROUTINE TPOUTU (CONTINUED) 
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CONT. ON PG 27 TPOUTJ 

PO 26 QF 22 

FIGURE 3.2.19. FUNCTIONAL FLOWCHART OF SUBROUTINE TPOUTJ (CONTINUED) 
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CGNT. GN PG 29 TPOUTJ 

PG 28 OF 22 — 

FIGURE 3.2.19. FUNCTIONAL FLOWCHART OF SUBROUTINE TPOUTJ (CONTINUED) 
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CONT. ON PG 31 


FIGURE 3.2.19. FUNCTIONAL FLOWCHART OF SUBROUTINE TPOUTJ (CONTINUED) 
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FORMAT C 1H+.84X. "NEVER ACTIVATED" ) 


CONT. ON PG 32 


TROUT J 


FIGURE 3.2.19. FUNCTIONAL FLOWCHART OF SUBROUTINE TPOUTJ (CONTINUED) 
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3,2.20 Subroutine: TREAD 


PURPOSE: 

METHOD: 

VARIABLES: 




Read the input timeline and control the creation of 
the event timeline. 

This routine reads an input timeline card, determines 
its type, and calls the correct routine to handle the 
type. 

The variables used in this subroutine are listed in the 
common blocks of the functional flowchart. Figure 3.2.20. 
See Appendix for definition of all variables. 
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FORMAT C2A1 . 17. 1 1 . A5. 2A2. A 1 . 2A5. 3X. 6A6 D 


ON PG 


TREAO 


FIGURE 3.2.20. FUNCTIONAL FLOWCHART OF SUBROUTINE TREAD 
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CQNT. ON PG 3 TREAD 

PG 2 QE 9 

FIGURE 3.2.20. FUNCTIONAL FLOWCHART OF SUBROUTINE TREAD (CONTINUED) 
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FIGURE 3 



CONT. ON PG 4 


TREAD 

P6 -3_— OE 9_ 


.2.20. FUNCTIONAL FLOWCHART OF SUBROUTINE TREAD (CONTINUED) 
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FIGURE 3.2.20. FUNCTIONAL FLOWCHART OF SUBROUTINE TREAD (CONTINUED) 
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CONI. ON PG 7 


TREAD 


FIGURE 3.2.20. FUNCTIONAL FLOWCHART OF SUBROUTINE TREAD (CONTINUED) 
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FIGURE 3.2.20. FUNCTIONAL FLOWCHART OF SUBROUTINE TREAD (CONTINUED) 
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WRITE C1U6.9G505 TYPE. NUMBER 


9050 

FORMAT C5X. 4 OHERROR ATTEMPTING TO REAO TIME LINE CARO . A2. [75 

3E 


I GO TO 10 



I GO TO 27Q > -4 
<1 


| WRITE C IUR5 ICOCC 1.25 

^5 


| 1PLC - ICOCC t. 2 5 


20Q_ 




^ A9 | 

i 



CONTINUE 

T 


| CALL RREADC IPLC.ODAT. 1000, ISTAT) 


f 


23 


CONT. ON PG 9 


TREAD 
P6. 8 . OF 


FIGURE 3.2.20. FUNCTIONAL FLOWCHART OF SUBROUTINE TREAD (CONTINUED) 


65 









210 




TREAD 


FIGURE 3.2.20. FUNCTIONAL FLOWCHART OF SUBROUTINE TREAD (CONTINUED) 
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3.3. PHASE II SUBROUTINES 
3.3.1 Subroutine: PHASE2 

PURPOSE: This routine controls the simulation of the Shuttle 

vehicle's electrical power system. 

METHOD: Using a user supplied time step this routine controls 

the following functions to simulate the vehicle from 
some input simulation start time to some input simulation 
abort time. 

1. Reads initialization data and solves the initialization 
calculations 

2. Determines the source I-V characteristics 

3. Solves the distribution system to find node 
voltages and branch currents 

4. Checks for constraint violations 

5. Provides the required outputs 

6. Repeats Steps 2 through 5 until the end of the 
simulation 

VARIABLES: The variables used in this subroutine are listed in the 
common blocks of the functional flowchart. Figure 3.3.1. 

See Appendix for definition of all variables. 
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SUBROUTINE PHASE2 
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c 

ALL COMMUNICATION IS TO TAKE PLACE THROUGH POP ELEMENTS I 
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STRAGl 
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ORIVER RELATED ELEMENTS 
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STRAG2 
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OC SYSTEM RELATED ELEMENTS 

c 

STRAG3 

ALL 

AC SYSTEM RELATED ELEMENTS 

c 

STRAG4 

ALL 

BATTERY RELATED ELEMENTS 

c 

STRAG5 

ALL 

FUEL CELL RELATED ELEMENTS 

c 

STRAG6 

ALL 

TR UNIT RELATED ELEMENTS 

c 

PROGRAM 

FLOW 



5 


INI TAL - READ OPTION CARD. FIXED DATA. INITIALIZE ALL SOURCES 
READ DOWN THE TIMELINE. DETERMINE INITIAL L0ADIN6 
CIRSQL - CIRCUIT SOLUTION DRIVER 

CALLS A. ACINVT - SOLVES AC INVERTERS APPLIES AS LOA 
TO DC CIRCUIT 

B. DCSGLV - SOLVES OC CIRCUIT 

C. TRAPLY - APPLIES TR UNIT AS LOAO TO AC CIRC 

D. ACSOLV - SOLVES AC CIRCUIT 

REDLIN - POWER MANAGEMENT 

INCRA - CONTROLS TIMESTEP IF REQUIRED REAOS TIMELINE INPUT 
CALLS A. TLREAD - READS CARD TIMELINE 
CHARGE - DETERMINES WHICH. IF ANY. 10 AH BATTERY IS CONNECTED 
TO THE BATTERY CHARGER AND UPDATES THAT BATTERYS STA 
QCLTMP - UPDATES BATTERY STATUS 

BATTIV - CALCULATES CORRECT BATTERY OPERATING IV CURVE 
FUCLTM - UPDATES FUEL CELL STATUS 

5 

FUCLIV - CALCULATES FUEL CELL OPERATING IV CURVE 

TRUNIT - CALCULATES TRANSFORMER-RECTIFIER UNIT OPERATING IV C 



CONT. ON PG 2 PHASE 2 

PG 1 OF a 


FIGURE 3.3.1 FUNCTIONAL FLOWCHART OF SUBROUTINE PHASE2 
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CONTINUE 


CALL CIRSGL 





CONT. ON PG 3 PHASE 2 

PG 2 OP 8 

FIGURE 3.3.1. FUNCTIONAL FLOWCHART OF SUBROUTINE PHASE2 (CONTINUED). 
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CONT. CN PG 4 


FIGURE 3.3.1. FUNCTIONAL FLOWCHART OF SUBRO 
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PHASE2 (CONTINUED) 






CGNT. ON PG 5 PHASE2 

PG 4 OF £ 

FIGURE 3.3.1. FUNCTIONAL FLOWCHART OF SUBROUTINE PHASE2 (CONTINUED) 
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FIGURE 3.3.1. FUNCTIONAL FLOWCHART OF SUBROUTINE PHASE2 (CONTINUED) 
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FIGURE 3.3.1. FUNCTIONAL FLOWCHART OF SUBROUTINE PHASE2 (CONTINUED) 
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FIGURE 3.3 


CONTINUE 
CALL TRUN1T 


EZJ 


I5Z 

"3“ 

£> 

AC 

, —l 

LZJ 


SZI 


At) 



CONTINUE 


1 GQ TO 1 Q 0l>t>g 



CONT. ON PG 8 

1 . FUNCTIONAL FLOWCHART OF SI 



PHASE2 (CONTINUED) 






[end [ 




PHASE 2 
P£_S FINAL 


FIGURE 3.3.1. 


FUNCTIONAL FLOWCHART OF SUBROUTINE PHASE2 (CONTINUED) 









3.3.2 Subroutine: ACINVT 


PURPOSE: 

METHOD: 


VARIABLES: 


**"*»*» Has 


To simulate the operation of the onboard dc-ac inverters 

The ac load and power factor for each inverter is 
calculated. These values are used to calculate the 
equivalent dc load and inverter efficiency. If the 
inverter is not carrying an ac load or the ac load it is 
carrying is less than the input dc no-load value, the 
inverter is set to dc no-load value. 

The variables used in this subroutine are listed in the 
common blocks of the functional flowchart. Figure 3.3.2 
See Appendix for definition of all variables. 
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CONT. ON PG 2 ACINVT 

PG 1 OF 6 

FIGURE 3.3.2. FUNCTIONAL FLOWCHART OF SUBROUTINE ACINVT 
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C0NT. ON PG 4 ACINVT 

PG 2 QE £ 

FIGURE 3.3.2, FUNCTIONAL FLOWCHART OF SUBROUTINE ACINVT (CONTINUED) 
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FIGURE 3.3.2. FUNCTIONAL FLOWCHART OF SUBROUTINE ACINVT (CONTINUED) 
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FIGURE 3.3.2. FUNCTIONAL FLOWCHART OF SUBROUTINE ACINVT (CONTINUED) 
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CALL INTERPC ACPFC I D , ACEFFC I ) . PFVC J ) . EFFT . 2 ) 


ACINt U=ACINLC I 3/ACEFFC I ) 


IFCMPRNT2.GT. OD 


WRITEC IU6. 901 OD [ . ACVAC I) . ACPFACC ID . AC INC ID 


CONT. ON PG 6 


ACINVT 


FIGURE 3.3.2. FUNCTIONAL FLOWCHART OF SUBROUTINE ACINVT (CONTINUED) 
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ACINVT 

P£_E FINAL 


FIGURE 3.3.2. FUNCTIONAL FLOWCHART OF SUBROUTINE ACINVT (CONTINUED) 
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3.3.3 Subroutine: BATTIV 

PURPOSE: To create the I-V curves used in the onboard battery 

simulations. 

METHOD: After determining the type of battery to be simulated, 

an interpolation is made into a set of current-voltage 
curves as a function of temperature to determine the 
battery I-V curves at its operating temperature. 

VARIABLES: The variables used in this subroutine are listed in the 

common blocks of the functional flowchart. Figure 3.3.3. 
See Appendix for definition of all variables. 
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C0NT. ON PG 2 


BATTIV 


FIGURE 3.3.3. FUNCTIONAL FLOWCHART OF SUBROUTINE BATTIV 
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FIGURE 3.3.3. 


CONT. ON PG 3 


BATTIV 
pg ? nr 


FUNCTIONAL FLOWCHART OF SUBROUTINE BATTIV (CONTINUED) 
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FIGURE 3.3.3. FUNCTIONAL FLOWCHART OF SUBROUTINE BATTIV (CONTINUED) 
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FIGURE 3.3.3. FUNCTIONAL FLOWCHART OF SUBROUTINE BATTIV (CONTINUED) 
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3.3.4 Subroutine: CHARGE 


PURPOSE: 

METHOD: 


VARIABLES: 


To simulate the operation of the battery charger. 

A charge curve is calculated based on depth-of-discharge 
and time since the last discharge. The curve is rate of 
charge versus time on charge. When the battery state-of- 
charge reaches 100 percent, the battery is removed from 
the charger. 

The variables used in this subroutine are listed in the 
common blocks of the functional flowchart, Figure 3.3.4. 
See Appendix for definition of all variables. 
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C SET ROUTING FLAGS SAVE DISTRIBUTION SWITCH POSITION 



CONT. ON PG 2 


CHARGE 
EG I OF 


FIGURE 3.3.4. FUNCTIONAL FLOWCHART OF SUBROUTINE CHARGE 
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CGNT. ON PG 3 CHARGE 

PS 2 P F 5 

FIGURE 3.3.4. FUNCTIONAL FLOWCHART OF SUBROUTINE CHARGE (CONTINUED) 
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CONT, ON PG 4 CHARGE 

pg 2 QE 5 

FIGURE 3.3.4. FUNCTIONAL FLOWCHART OF SUBROUTINE CHARGE (CONTINUED) 
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CGNT. ON PG 5 * 


CHARGE 
PG 4 OF 


FIGURE 3.3.4. FUNCTIONAL FLOWCHART OF SUBROUTINE CHARGE (CONTINUED) 
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FIGURE 3.3.4. FUNCTIONAL FLOWCHART OF SUBROUTINE CHARGE (CONTINUED) 
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3.3.5 Subroutine: DCSOLV 


PURPOSE: 

METHOD: 

VARIABLES: 


To determine the operating characteristics of the 
EPDC. 

Utilizing a user provided distribution circuit this 
routine calculates the nodal equations of the circuit 
and solves these equations by means of the Gauss-Jordan 
method. This process is done iteratively until the 
change-in node voltages are less than some input value. 

The variables used in this subroutine are listed in the 
common blocks of the functional flowchart. Figure 3.3.5. 
See Appendix for definition of all variables. 
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subroutine dcsolv 

INCLUDE STRAG1 
INCLUDE STRAG2 
DIMENSION AMC 12) 
DIMENSION AVU2) 
DIMENSION BVC 12) 
DIMENSION II C 12) 
DIMENSION STAMC12) 


DIMENSION VO(l 2) 
DIMENSION VSCI2) 
DIMENSION VS0C12) 
DIMENSION BC30) 
DIMENSION QC30.30) 
DIMENSION WC30) 
DIMENSION YC30.30) 
DATA CST 1 / l.E-23/ 


DATA CST2 /-l.EIO/ 
DATA CST3 / l.E-10/ 

Data cst4 / i.e-?/ 

DATA ES /12*28./ 
DATA NOR /{/ 

DATA MPRINT /0/ 

DATA NITER /30/ 
ITER=0 


COMPUTE OPEN CIRCUIT VOLTAGES 


< DO 100 iVl .NSCS > - 


> 2 


JFt NSOURC I . 1 ) . GT. NCRT ) 

' T 

L GO TO lOQ 


CALL INTERPt 0. . VOt I ). SCC 1.2.1). SCC 1.1.1). NSCC I ) ) 


r 




CONT. ON PG 2 


DCSOLV ’ 

ee_ 1 DE- 




FIGURE 3.3.5. FUNCTIONAL FLOWCHART OF SUBROUTINE DCSOLV 
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CONT. ON PG 3 DC SOL V 

PG 2 OP 22- 

FIGURE 3.3,5. FUNCTIONAL FLOWCHART OF SUBROUTINE DCSOLV (CONTINUED) 
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N2=ICRCT( 18.2) 

VLCn=SCIB}*CVCNl)-VCN2)-CRCTC IB. 1 )-CBC IB)*CRC IB)+CRCTC IB.5))) 
ZX= 1 . E+20 
ZY=I .E+20 



FIGURE 3.3.5. FUNCTIONAL FLOWCHART OF SUBROUTINE DCSOLV (CONTINUED) 
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FIGURE 3.3.5. FUNCTIONAL FLOWCHART OF SUBROUTINE DCSOLV (CONTINUED) 
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C FIND THE BOUNDING INTERVAL ON THE SOURCE CURVE 


IB=NS0UR( I . n 


CONT. ON PG G 


DC SOL V 

PG._5 QE 22- 


FIGURE 3.3.5. FUNCTIONAL FLOWCHART OF SUBROUTINE DCSOLV (CONTINUED) 
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FIGURE 3.3.5. FUNCTIONAL FLOWCHART OF SUBROUTINE DCSOLV (CONTINUED) 
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CONTINUE 



/ 

FIGURE 3.3.5. FUNCTIONAL FLOWCHART OF SUBROUTINE DCSOLV (CONTINUED) 
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CONT. ON PG 10 


FIGURE 3.3.5. FUNCTIONAL FLOWCHART OF 



DCSOLV (CONTINUED) 








CONT. ON PG ! 1 


OCSGUV 


FIGURE 3.3.5. FUNCTIONAL FLOWCHART OF SUBROUTINE DCSOLV (CONTINUED) 
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CONT. ON PG 12 OCSOLV 

PG 11 0 £ 21 

FIGURE 3.3.5. FUNCTIONAL FLOWCHART OF SUBROUTINE DCSOLV (CONTINUED) 
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FIGURE 3.3.5. FUNCTIONAL FLOWCHART OF SUBROUTINE DCSOLV (CONTINUED) 
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CGNT. ON PG 14 


DCSOLV 


FIGURE 3.3.5. FUNCTIONAL FLOWCHART OF SUBROUTINE DCSOLV (CONTINUED) 









FIGURE 3.3.5. FUNCTIONAL FLOWCHART OF SUBROUTINE DCSOLV (CONTINUED) 


211 

















1FCMPRINT.LT. 1 D 


FIGURE 3. 



I GO TO 156(T 

rd- 


~ < DO 1550 J=1 .N'NODE~> 
_s 


WRtTEC I U6. 50201 CQU.JT. I=1.NN0DE1 
VRITEC 1UG. 5021 ) WCJ) 


. 5Q 2 0 $ 

j FQRMATC5X. •'Q”»4X. 10F10.4) 


502 T ± 


FQRMATC” II4X. V" . 4X. F 1 0. A // 1 


1 55.0 A 


CONTINUE 


1 56 0 ± 


CONTINUE 


ZERO THE REFERENCE NOOE 
- 

<00 1600 1 . NNODF > - 

'if 


- - >16 


OCNOR. I )=0. 
QCl.N0R5=0. 


I 


CONT. ON PG 16 


OCSOLV 

ES..15.-QF 22 - 


J.5. FUNCTIONAL FLOWCHART OF SUBROUTINE DCSOLV (CONTINUED) 
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FIGURE 3.3.5. FUNCTIONAL FLOWCHART OF SUBROUTINE DCSOLV (CONTINUED) 
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FORMAT C5X. “NODE VALUES LAST I TERAT IQN'VC 5X. 1 OF L 0. 4) ) | 



CONT. ON PG 18 


DCSOLV 


FIGURE 3.3.5. FUNCTIONAL FLOWCHART OF SUBROUTINE DCSOLV (CONTINUED) 
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FIGURE 3.3.5. FUNCTIONAL FLOWCHART OF SUBROUTINE DCSOLV (CONTINUED) 
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FIGURE 
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' " <00 i fl?0 T=1 , U££sO > 
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IB-NSOURC I * t ) 

WRITEC I US* 5040) I . CSC I). ESC I ). VSC I ).St 18) 


5040 6 

F0RMATC5X* " I CS ES VS S ,f . 1 2. 3F 1 0.4. 15) 


CONTINUE 


VRITEC IU6.5041) TDEL. DELTA 


5041 6 

FORMAT C /5X • ’■ TDEL DELTA" . 5X* 2F 15. 5// 3 

I WRITEC IUS. 50503 C I . VLC I } . CLC I) ► I=i.NLDS) 
* ■ 1 " * ..... . ... 


❖ 

CGNT. ON PG 20 DCSOLV 

PG 19 OF 27 


.3.5, FUNCTIONAL FLOWCHART OF SUBROUTINE DCSOLV (CONTINUED) 
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TEST FOR REVERSE BIAS ON THE BRANCHES CONTAINING OlODES 



FIGURE 3.3.5. FUNCTIONAL FLOWCHART OF SUBROUTINE DCSOLV (CONTINUED) 
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IFCCBt I }. GT. -CST4} 




CONT. ON PG 22 

FIGURE 3.3.5. FUNCTIONAL FLOWCHART OF SUBROU 
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DCSOLV (CONTINUED) 











WRITECIUG.5055) I.NI.N2.V(N1).V(N2). CRCT C 1.2). 
CRCT(I.3).CRCTC 1.5) 



CQNT. ON P6 23' 


DCSOLV 


FIGURE 3.3.5. FUNCTIONAL FLOWCHART OF SUBROUTINE DCSOLV (CONTINUED) 

7 


219 


















CGNT. ON PG 24 


OCSOLV 

PG 23 ‘OP 27 


FIGURE 3.3.5. FUNCTIONAL FLOWCHART OF SUBROUTINE DCSOLV (CONTINUED) 
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CONT. ON PG 25 


OCSOLV 


FIGURE 3.3.5. FUNCTIONAL FLOWCHART OF SUBROUTINE DCSOLV (CONTINUED) 
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CONTINUE 


CONT. ON PG 26 


DCSOLV 


FIGURE 3.3.5. FUNCTIONAL FLOWCHART OF SUBROUTINE DCSOLV (CONTINUED) 
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FIGURE 3.3.5. FUNCTIONAL FLOWCHART OF SUBROUTINE DCSOLV (CONTINUED) 
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3.3.6 Subroutine: FUCLIV 


PURPOSE: 

METHOD: 

VARIABLES: 


To create the I-V curves used in simulating the onboard 
fuel cells 

An interpolation is made into a set of current-vol tage 
curves as a function of temperature to determine the fuel 
cell I-V curve at its operating temperature. The curve 
is further degraded by the parasitic load to be carried 
by the fuel cel 1 . 

The variables used in this subroutine are listed in 

the common blocks of the functional flowchart, Figure 3.3.6. 

See Appendix for definition of all variables. 
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FUCLIV 
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FIGURE 3.3.6. FUNCTIONAL FLOWCHART OF SUBROUTINE FUCLIV 
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j FCCVCK. 2) = 0. § 
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INUE 
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H — | ao 
20 2 _ 


CONTINUE 

T 


{ CALL INTRPJCFCCVC 1. 1) . FCTC J ) . FCTN. 7. FCTAC 1 . 2) . NFCTa:) 


=3 


' < 00 50 K=1 .NFCTA > 


FCCVCK. 2)= 

FCTACK. 1 ) 

- - 50 ! 
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| CONTINUE 




CONT. ON PS 3 


FUCLIV 
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FIGURE 3,3.6. FUNCTIONAL FLOWCHART OF SUBROUTINE FUCLIV (CONTINUED) 
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FIGURE 3.3.6. FUNCTIONAL FLOWCHART OF SUBROUTINE FUCLIV (CONTINUED) 
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FIGURE 3.3.6. FUNCTIONAL FLOWCHART OF SUBROUTINE FUCLIV (CONTINUED) 
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3.3.7 Subroutine: FUCLTM 


PURPOSE: To determine the change in fuel cell temperature and 

certain gross cyrogenic quantities. 

METHOD: Based on the current fuel cell operating temperature 

and steady state temperatures versus current curves an 
ideal operating current is determined. The difference 
between the ideal current and the operating current is 
used to determine the change in fuel cell temperature. 

Based on input purge times and rates, and input usage 
rates, the amount of oxygen and hydrogen used and the 
quantity of water produced is calculated. 

VARIABLES: The variables used in this subroutine are listed in the 

conrnon blocks of the functional flowchart. Figure 3.3.7. 
See Appendix for definition of all variables. 
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SUBROUTINE FUCLTMC J 3 
INCLUDE STRAG1 
INCLUDE STRAG5 


THIS SUBROUTINE PERFORMS THE FOLLOWING FUEL CELL CALCULATIONS 
HYOROGEN AND OXYGEN DEPLETION 

WATER PRODUCTION - 

HYDROGEN AND OXYGEN PURGES 
CRYOGENIC DEPLETION 
INPUT 

DTM TINE STEP 

FCCP OPERATING CURRENT QF FUEL CELL 


HYOROGEN PURGE RATE 
TINE REQUIREO FOR HYDROGEN PURGE 
HYOROGEN USAGE RATE 
FUEL CELL NUMBER 
OXYGEN PURGE RATE 
TINE REOUIREO FOR OXYGEN PURGE 
OXYGEN USAGE RATE 
INTERVAL BETWEEN HYDROGEN PURGES 


PIO INTERVAL BETWEEN OXYGEN PURGES 

TMAXFC MAX ALLOWABLE TIME FOR ONE PASS THRU 
TPH TIME OF LAST HYDROGEN PURGE 

TPO TIME OF LAST OXYGEN PURGE 

WPR WATER PRODUCTION RATE 

INPUT/OUTPUT 

H20T TOTAL AVAILABLE WATER 

H2T TOTAL AVAILABLE HYDROGEN 


TOTAL AVAILABLE OXYGEN 
OUTPUT 

TIME FUEL CELL HAS BEEN OPERATING 
HYOROGEN DEPLETED THIS T [MESTEP 
OXYGEN DEPLETED THIS TIMESTEP 
WATER PRODUCED THIS TIMESTEP 


COMPUTE INTERNAL TIME VARIABLES 

KMAX^CDTM/TMAXFU+U 

QTFC=OTM/KMAX 


ON PG 


FUCLTM 


[GURE 3.3.7. FUNCTIONAL FLOWCHART OF SUBROUTINE FUCLTM 
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FIGURE 3.3.7. FUNCTIONAL FLOWCHART OF SUBROUTINE FUCLTM (CONTINUED) 
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CONT. ON PG 4 
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FIGURE 3.3.7. FUNCTIONAL FLOWCHART OF SUBROUTINE FUCLTM (CONTINUED) 
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H2T=H2T-H2 
H2GT =H20T +H20 
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FIGURE 3.3.7. FUNCTIONAL FLOWCHART OF SUBROUTINE FUCLTM (CONTINUED) 
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FIGURE 3.3.7. FUNCTIONAL FLOWCHART OF SUBROUTINE FUCLTM (CONTINUED) 
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FIGURE 3.3.7. FUNCTIONAL FLOWCHART OF SUBROUTINE FUCLTM (CONTINUED) 
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3.3.8 Subroutine: INCRA 

PURPOSE: To control the sequential time dependent operation of 

Phase II. 

METHOD: The present time {MET or GET) is incremented by the 

input time step. Both the input card timeline and the 
interface tape are checked to see if either or both 
should be read. If both are to be read, the interface 
tape is read first, then the card timeline. 

VARIABLES: The variables used in this subroutine are listed in the 

common blocks of the functional flowchart. Figure 3.3.8. 
See Appendix for definition of all variables. 
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COMMON /TQTPWR/ PWRTQT 

COMMON /PHAS 1 / I DUM{ 3855 ).T ABORT 

DATA TLOAD / 100*0./ 

DATA SUBSYS / 100*0./ 

DATA NIF /0/ 

DATA IARR1 /52. 61 . 62. 63. 53. 64. 65. 66. 54. 67. 68. 69/ 

6 

C INCREMENT THE PRESENT TIME 

. ' 
SMET=MET 
MET=MET+TOELTA 

, 3E 

C CHECK FOR TIMELINE 

r . 

YTIME=0. 


IFCTREADC+. 00 1 . LT. TREADT ) 


CONT. ON PG 2 


INCRA 


FIGURE 3.3.8. FUNCTIONAL FLOWCHART OF SUBROUTINE INCRA 
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FIGURE 3.3.8. FUNCTIONAL FLOWCHART OF SUBROUTINE INCRA (CONTINUED) 
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CONT. ON PG 4 INCRA 

PG 3 OF 8 



FIGURE 3.3.8. FUNCTIONAL FLOWCHART OF SUBROUTINE INCRA (CONTINUED) 
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FIGURE 3.3.8. FUNCTIONAL FLOWCHART OF SUBROUTINE INCRA (CONTINUED) 
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FIGURE 3.3.8. FUNCTIONAL FLOWCHART OF SUBROUTINE INCRA (CONTINUED) 
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A9 


IF C IFLG05.LT. n 



j CONTINUE I 

CALL MCGMP2CSMET.TDUM.2.TLQA0.SU8SYS. li 


IFCNIF+QQASGC IU83.EG. OD 


CGNT. ON PG 7 


INCRA 


FIGURE 3.3.8. FUNCTIONAL FLOWCHART OF SUBROUTINE INCRA (CONTINUED) 
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CALL CTAP2CMET. TDELTA. PCDCT C H . ICN. IU8* 

PFB.S) 



FIGURE 3*3.8, FUNCTIONAL FLOWCHART OF SUBROUTINE INCRA (CONTINUED) 
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3.3.9 Subroutine: INITAL 

PURPOSE: To prepare Phase2 for execution. 

METHOD: All initialization procedures are followed and all 

necessary initial quantities are calculated. 

VARIABLES: The variables used in this subroutine are listed in the 

common blocks of the functional flowchart. Figure 3.3.9. 
See Appendix for definition of all variables. 


preceding page blank trnjsm 
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C0NT. ON PG 2 INITAl 

PG._1_.QF__ . .10 

FIGURE 3.3.9. FUNCTIONAL FLOWCHART OF SUBROUTINE INITAL 
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CONI. ON PC 3 


FIGURE 3.3.9. FUNCTIONAL FLOWCHART OF SUBROUTINE INITAL (CONTINUED) 
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CQNT. ON PG 4 


FIGURE 3.3.9. FUNCTIONAL FLOWCHART OF SUBROUTINE INITAL (CONTINUED) 
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POOR QUALITY 


252 















FIGURE 3.3.9. FUNCTIONAL FLOWCHART OF SUBROUTINE INITAL (CONTINUED) 
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FIGURE 3.3.9. FUNCTIONAL FLOWCHART OF SUBROUTINE INITAL (CONTINUED) 
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50 & 

CONTINUE 

V 

READC IU5.9020,ENO=SO: A 
WRITECIU6 .9020) A 
WRITECIUI 1 . 9020) A 

^ ' 

I GO TO 50 > 



FIGURE 3.3.9. FUNCTIONAL FLOWCHART OF SUBROUTINE INITAL (CONTINUED) 
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FIGURE 3.3.9. FUNCTIONAL FLOWCHART OF SUBROUTINE INITAL (CONTINUED) 


ORIGINAL PAGE IS 

of poor quality: 
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FIGURE 3.3.9. FUNCTIONAL FLOWCHART OF SUBROUTINE INITAL (CONTINUED) 
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FIGURE 3.3.9. FUNCTIONAL FLOWCHART OF SUBROUTINE INITAL (CONTINUED) 


OKieiNAIJ PAGE IS 
OF POOR QUALITY 









FIGURE 3. 



3QM 

FORMAT ( 1H1 ) 


9020 v 

FORMAT ( 8 0A13 


,9Q3Q i 

FORMATCFIO.O) 


END 


INITAL 

PG 10 FINAL 


.9. FUNCTIONAL FLOWCHART OF SUBROUTINE INITAL (CONTINUED) 
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3.3.10 Subroutine: QCLTMP 


PURPOSE: To update battery status 

METHOD: Based upon time elapsed, operating temperature, and 

operating current and voltage the following quantities 
are calculated: 

1. Charge and discharge efficiency 

2. Heat produced 

3. Change in capacity 

4. Temperature 

5. New ampere-hour status 

6. State-of-charge 

VARIABLES: The variables used in this subroutine are listed in the 
common blocks of the functional flowchart. Figure 3.3.10. 
See Appendix for definition of all variables. 
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n n no no o n 


subroutine qcltmpci. IACT) 

INCLUDE STRAG1 
INCLUOE STRAG2 
INCLUDE STRAG4 

e: ~ : 

PRIMARY BATTERY SUBROUTINE CALCULATES A-H STATUS. SOC. OEGRAT I ON 
TURNS ON HEATERS. CALCULATES HEAT GENERATED. AND TAKES OEPLETEO 
BATTERIES OUT OF LINE 

INPUT 

CS OPERATING CURRENT OF 8ATTERY 

CSUBD BATTERY DESIGN CAPACITY 

DT TIMESTEP 

EFF DISCHARGING BATTERY EFFICIENCY 


EFFAVR CHARGING BATTER EFFICIENCY 

PUP2.P3.P4 CONSTANTS FOR HEAT GENERATION EQUATION 
OHI HEATER OUTPUT CAPABILITY (WATTS) 

TB BATTERY TEMPERATURE 

VP OPERATING VOLTAGE OF BATTERY 

TRFCCJ.H TEMP VS TEMP RATE OF INCREASE 
1=1 TEMP RATE FOR CURRENT CDX1 
1=2 TEMP RATE FOR CURRENT CDX2 


c 


1=3 TEMP RATE FOR CURRENT CDX3 

c 


I -A TEMP IN DEGREES F 

c 


' </= POINTS 

c 


OUTPUT 

c 

AH 

AMP HOUR STATUS OF BATTERY 

c 

CC 

AMP HOUR CAPACITY OF BATTERY 

c 

00 

HEAT PRODUCED CWATTS) 

c 

QH 

HEATER OUTPUT C WATTS) 



C 

SOC 

PRESENT 

STATE-OF-CHARGE OF 

THE BATTERY 

c 

TD 

TIME ON 

LAST MAJOR BATTERY 

DISCHARGE 


DIMENSION TRFCC8.4) 

OATA C0X1.CDX2.C0X3 /I5. . 25. 8. 36. 7/ 

■■■■ — g 

DATA TRFC /7..G. 6. 5. 4. 3. 8. 1.5. -I. 2. -5. .-9. 4. 

9.2*9. .8. 1*6. 8*5. .2.7. 0, .—3*8. 
12.6. 12.5. 12.. I 1.2. 9. 8. 8. .5.9.3.. 
30. .40. .60. .80. . 100. . 120. . 140. . ISO. 


CONT. ON PG 2 


GCLTMP 


FIGURE 3.3.10. FUNCTIONAL FLOWCHART OF SUBROUTINE QCLTMP 
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C0NT. ON PG 3 QCLTMP 

PG 7 OF B 

FIGURE 3.3.10. FUNCTIONAL FLOWCHART OF SUBROUTINE QCLTMP (CONTINUED) 
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C CALC DESCENT 8ATT TEMP 



CONT. ON PG 4 QCLTMP 

P5 3—SE 6 

FIGURE 3.3.10. FUNCTIONAL FLOWCHART OF SUBROUTINE QCLTMP (CONTINUED) 
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FIGURE 3.3.10. 


CONT. ON PG 5 


QCLTNP 

PG 4 OF 6 


FUNCTIONAL FLOWCHART OF SUBROUTINE QCLTMP (CONTINUED) 
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C0NT. ON PG 6 QCLTMP 

P5 5 OF 6 

FIGURE 3.3.10. FUNCTIONAL FLOWCHART OF SUBROUTINE QCLTMP (CONTINUED) 
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END 


ORIGINAL] PA0B13 

^ POOR QUALITY 


OCLTHP 
PG 6 FINAI 


FIGURE 3.3.10. FUNCTIONAL FLOWCHART OF SUBROUTINE QCLTMP (CONTINUED) 
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3.3.11 Subroutine: REDLIN 


PURPOSE: To detect and identify EPS values which violate established 

limits 

METHOD: The subroutine consists of logic for tests on six individual 

EPS parameters. The parameters are: 

1. Node Voltage - tested for violation of under voltage limit 

2. Inverter Volt-amps - checked for overload 

3. Branch Current - tested for current limit 

4. Fuel Cell Power 

a. Peak Power - tested for power level exceeding peak 

power limit 

b. Continuous Power - tested for power level exceeding 

Continuous power limit 

c. Minimum Power - tested for power level under minimum 

power level 

5. Cryogen Level 

a. O 2 - checked for depletion 

b. H 2 - checked for depletion 

6. Battery State-of-Charge - checked for depletion 

The tests are made by comparing an EPS parameter value to its 
appropriate limit from the fixed data tape. If the parameter 
value violates the limit, a message is constructed which 
contains the following data: 

1. Time of Violation 

2. Type of Violation 

3. Value of Limit 

4. Parameter Value 

VARIABLES: The variables used in this subroutine are listed in the common 

blocks of the functional flowchart, Figure 3.3.11. See Appendi 
for definition of all variables. 
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FIGURE 3.3. 


SUBROUTINE REDlINCIUJ 
INCLUDE STRAG1 
INCLUDE STRAG2 
INCLUDE STRAG3 
INCLUDE STRAG4 
INCLUDE STRAG5 
INCLUDE STRAG6 
INCLUDE STRAG7 
•6 


DIMENSION KFLGC3.5J 
DIMENSION SFCPC3.53 
DIMENSION AVPWRC3.5) 
DIMENSION ATIMEC3.5J 
DIMENSION VTIMC3.5D 
DIMENSION IS0CC6J 
DATA KFLG /15*0 / 
DATA ISOC / 6*0 / 


< 00 100 1=1 , NNODE > 



I GO TO 1 00 


ITYPEV=1 


, v 

| WRITEUU3 MET. ITYPEV. I.UVCn.VCI) 


Iflfl 4 

CONTINUE 

i 


< no 700 r=I .NfNV > - - - 

V 


> 2 


* 

com. ON PG 2 REDLIN 

P6 0£ 9. 


11. FUNCTIONAL FLOWCHART OF SUBROUTINE REDLIN 
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FIGURE 3.3.11. FUNCTIONAL FLOWCHART OF SUBROUTINE REDLIN (CONTINUED) 
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FIGURE 3.3.11. FUNCTIONAL FLOWCHART OF SUBROUTINE REDLIN (CONTINUED) 
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FIGURE 3.3. 




v 

CONT. ON PG 5 REDLIN 

PG 4 OF 9 


11. FUNCTIONAL FLOWCHART OF SUBROUTINE REDLIN (CONTINUED) 
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CCNT. ON PG 6 


REDLIN 
PG 5 OF 


FIGURE 3.3.11. FUNCTIONAL FLOWCHART OF SUBROUTINE REDLIN (CONTINUED) 











CGNT. ON PG 7 


REOLIN 


FIGURE 3.3.11. FUNCTIONAL FLOWCHART OF SUBROUTINE REDLIN (CONTINUED) 
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CONT. ON PG 8 


REOLIN 


FIGURE 3.3.11. FUNCTIONAL FLOWCHART OF SUBROUTINE REDLIN (CONTINUED) 
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CONT. ON PG 9 


REDLIN 


FIGURE 3.3.11. FUNCTIONAL FLOW CHART OF SUBROUTINE REDLIN (CONTINUED) 
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REDLIN 

PG 3 FINAL 


FIGURE 3.3.11. FUNCTIONAL FLOWCHART OF SUBROUTINE REDLIN (CONTINUED) 
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3.4 ANALYSIS SUBROUTINES 
3.4.1 Subroutine: COMUSE 


PURPOSE: To provide a component analysis 

METHOD: Using the component dictionary, compacted dictionary, and 

interface tape the following quantities are determined 
for: 

1. Mission Phase: 

Subsystem Analysis 

Average kilowatts (KW) 

Percent of total KW 
Maximum KW 
Time of maximum 

Phase Analysis 

Total kilowatt hours (KWH) required 

Accumulated KWH 

Average KW for mission phase 

Maximum KW 

Time of maximum (KW) 

2. Component 

Average use factor 
Total on time 

Component energy required (WH) 

Percent of total mission energy required 

3. Subsystem 

By component 

Average use factor 
Total on time 

Component energy required (WH) 

Percent of subsystem energy required 

VARIABLES: The variables used in this subroutine are listed in the 

common blocks of the functional flowchart. Figure 3.4.1. 
See Appendix for definition of all variables. 
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CQNT. ON PG 2 COMUSE 

PG 1 OP L2 

FIGURE 3.4.1. FUNCTIONAL FLOWCHART OF SUBROUTINE COMUSE 


• original page is 

OF POOH QUALITY 
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CONT. ON PG 3 COMUSE 

P 5 2 Q.F 12 

FIGURE 3.4.1. FUNCTIONAL FLOWCHART OF SUBROUTINE COMUSE (CONTINUED) 
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FIGURE 3.4.1. FUNCTIONAL FLOWCHART OF SUBROUTINE COMUSE (CONTINUED) 
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CGNT. ON PG 5 COMUSE 

PS 4 OF L2_ 

FIGURE 3.4. 1,. FUNCTIONAL FLOWCHART OF SUBROUTINE COMUSE (CONTINUED) 
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| CALL liCOMUSC T2. 2. PLOAD. LNUM< PFAC ) 




CONT. ON PG 6 COMUSE 

P3 5 OF 12- 

FIGURE 3.4.1 . FUNCTIONAL FLOWCHART OF SUBROUTINE COMUSE (CONTINUED) 
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FIGURE 3 . 4 . 1 . FUNCTIONAL FLOWCHART OF SUBROUTINE COMUSE (CONTINUED) 












FIGURE 3.4.1. FUNCTIONAL FLOWCHART OF SUBROUTINE COMUSE (CONTINUED) 
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< 00 80 1 = 1. tCN > - - - t> 9 
V 

CONT. ON PG 9 COMUSE 

P5 8 OP L2 

FIGURE 3.4.1. FUNCTIONAL FLOWCHART OF SUBROUTINE COMUSE (CONTINUED) 


DRIGIN4BvPAGB is 
of pgoh quality 


287 





















288 



COMUSE (CONTINUED) 










CONT. ON PG 11 COM USE 

P3 10 OF 12 


FIGURE 3.4.1. FUNCTIONAL FLOWCHART OF SUBROUTINE COMUSE (CONTINUED) 
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CONT. ON PG 12 COMUSE 

P3 n OF L2 


FIGURE 3.4.1. FUNCTIONAL FLOWCHART OF SUBROUTINE COMUSE (CONTINUED) 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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COMUSE 


FIGURE 3.4.1. FUNCTIONAL FLOWCHART OF SUBROUTINE COMUSE (CONTINUED) 
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APPENDIX 



APPENDIX 


PKECEDItiQ 





SEPS DATA VARIABLES 


PDP COMMON VARIABLE 

ELEMENT BLOCK NAME l_ J_ K 

STRAG1 CONTRL DT 

IFLGOl 

IFLG02 

IFLG03 

IFLG04 

IFLG05 

IFLG06 

IFLG07 

IFLG08 

IFLG09 

IFLGIO 

IFLG11 

IFLG12 

IFLG13 

IFLG14 

IFLG15 

IFLG16 

IFLG17 


DEFINITION 

Actual time interval between 
successive time steps in decimal 
hours 

Option Flag 

> 0 - Execute Phase I 

= 0 or blank - Do not execute 
Phase I 

> 0 - Execute Phase II 

= 0 or blank - Do not execute 
Phase II 

> 0 - Plot on Printer 1 

= 0 or blank - Do not plot on 
Printer 1 

Not Used 

> 0 - Use 3 point load data 

s 0 or blank - Do not use 3 point 
load data 

Not Used 
Not Used 
Not Used 
Not Used 

> 0 - Print each Phase II timepoint 
= 0 or blank - Do not print each 

Phase II timepoint 

Not Used 
Not Used 
Not Used 
Not Used 
Not Used 
Not Used 
Not Used 


A-l 



SEPS DATA VARIABLES 


POP COMMON VARIABLE 

ELEMENT BLOCK NAME I_ J 

STRAG1 CONTRL IFLG18 

IFLG19 

IFLG20 

IFLG21 

IFLG22 

IFLG23 

IFLG24 

IFLG25 

IFLG26 

IFLG27 

IFLG28 

IFLG29 

IFLG30 

IFLG31 

IFLG32 

IFLG33 

IFLG34 

IFLG35 


K_ DEFINITION 

Not Used 
Not Used 

> 0 - Execute Phase II COMUSE 
= 0 or blank - Do not execute 

COMUSE 

Not Used 
Not Used 
Not Used 
Not Used 
Not Used 
Not Used 
Not Used 
Not Used 

> 0 - Suppress analysis part 1 
= 0 or blank - Do not suppress 

analysis part 1 

> 0 - Execute Phase I COMUSE 
= 0 or blank - Do not execute 

Phase I COMUSE 

> 0 - Print out input components 
= 0 or blank - No action taken 

> 0 - Print out input procedures 
= 0 or blank - No action taken 

> 0 - Print out input activities 
= 0 or blank - No action taken 

> 0 - Print out input timeline 
= 0 or blank - No action taken 

> 0 - Suppress COMUSE component 

analysis 

- 0 or blank - Do not suppress 
COMUSE component 
analysis 


A-2 



SEPS DATA VARIABLES 


POP COMMON VARIABLE 

ELEMENT BLOCK NAME I_ J_ K 

STRAG1 CONTRL IFLG36 

IFLG37 

IFLG38 

IFLG39 

IFLG40 
I ABORT 

IFILE 20 


DEFINITION 

> 0 - At the end of each mission 

phase provide a listing of 

all active components 

= 0 or blank - No action taken 

> 0 - Mission phase analysis 

- 0 or blank - No mission phase 
analysis 

> 0 - Suppress cycled component 

listing 

s 0 or blank - Do not suppress 
cycled component 
1 i sti ng 

> 0 - Suppress print of compacted 

dictionary 

5 0 or blank - Do not suppress 

compacted dictionary 
printout 

> 0 - Suppress subsystem analysis 

at each time point 

* 0 or blank - Do not suppress 

subsystem analysis 

Set > 0 to abort simulation 


1=1 File containing components 

1=2 File containing procedure 

1=3 File containing activities 

1=4 File containing timelines 

1=5 File containing fixed data 

1=6 Not Used 

1=7 Not Used 

1=8 File containing Phase I interfaces 

1=9 File containing Phase I plot 

1=10 File containing Phase II interface 


A- 3 



SEPS DATA VARIABLES 


PDP 

COMMON 

VARIABLE 


ELEMENT 

BLOCK 

NAME I J K 

DEFINITION 

STRAG1 

CONTRL 

IFILE 1=11 

File containing Phase II plot 



1=12 

File containing compacted 
dictionary 

V 



1=13 

File containing circuit 
definition 



1=14 

File containing constraint 
viol ations 



1=15 

File containing subsystem 
names 



1=16 

File containing mission phase 



1=17 

Not Used 



1=18 

Not Used 



1=19 

Not Used 



1=20 

Not Used 


ISOLVC Flag to request a circuit solution 

at a particular card timeline point. 
Set = 1 to request solution 

IUNIT 20 

1=1 Unit containing components 

1=2 Unit containing procedures 

1=3 Unit containing activities 

1=4 Unit containing timeline 

1=5 Unit containing fixed data 

1=6 Not Used 

1=7 Not Used 

1=8 Unit containing Phase I 
interface 

1=9 Unit containing Phase I plot 


A- 4 



SEPS DATA VARIABLES 


PDP COMMON VARIABLE 

ELEMENT BLOCK NAME 


STRA61 CONTRL IUNIT 


J PRINT 

NPRT 

MET 
TDELTA 
TREADC 
TREADT 
DEBUG MPRNT 
MPRNT1 
MPRNT2 


J K DEFINITION 


1=10 

Unit containing 
interface 

Phase II 

1=11 

Unit containing 
plot 

Phase II 

1=12 

Unit containing 
dictionary 

compacted 

1=13 

Unit containing 
definition 

circuit 

1=14 

Unit containing 
violations 

constraint 

1=15 

Unit containing 
names 

subsystem 

1=16 

Unit containing 
phases 

mission 

1=17 

Not Used 


1=18 

Not Used 


1=19 

Not Used 


1=20 

/ 

Not Used 



Flag to request initialization 
data to be printed out. Value >1 

Formatted printout interval as a 
multiple of TDELTA 

Simulation start time 

Maximum simulator time increment 

Next time to read card input 

Next time to read tape input 

Debug print control 

Debug print control 

Debug print control 


A- 5 



SEPS DATA VARIABLES 


PDP 

ELEMENT 

COMMON 

BLOCK 

VARIABLE 

NAME 

I J K DEFINITION 

ST RAG! 

DEBUG 

MPRNT3 

Not Used 

- 


MPRNT4 

Debug print control 



MPRNT5 

Not Used 



MPRNT6 

Not Used 


CMPCNT 

CVAL 

25 

Component power value to be used in 
Phase II calculations in lieu of 
the component loads in TPLOAD. 

NOTE: LOC 1-5 for F/C 1-5 noncyclic 

loads and LOC 6-10 for F/C 1-5 
cyclic loads 

CMPCNT 

NCNT 

25 

Component numbers associated with 
the component loads in CVAL 


NCNTC 


Number of loads defined in CVAL 


NCTP 

25 

For each LOAD in CVAL, defines if 
the load is constant power (>0) or 
constant resistance (=0) 

UNITS 

IU5 


Internal unit designation 


IU6 


Internal unit designation 


IU7 


Internal unit designation 


IU8 


Internal unit designation 


IU9 


Internal unit designation 


IU10 


Internal unit designation 


IU11 


Internal unit designation 

CONSTR 

ACVA 

9 

AC volt-ampere load on inverter 


CAPINV 

9 

Inverter (I) maximum overload limit 


FCLIM 

3 

Fuel cell power limits, 1 = peak, 
2 = average, 3 = minimum 


H2U 


Unusable hydrogen quantity 


02U 


Unusable oxygen quantity 


A- 6 



SEPS DATA VARIABLES 


PDP 

ELEMENT 

COMMON 

BLOCK 

VARIABLE 

NAME 

I_ J_ 

K DEFINITION 

STRAG1 

CONSTR 

SOCUL 


Minimum battery SOC limit 



TIMV 

3 

Length of time that fuel cell power 
limits apply, 1 = peak, 2 = average, 
3 = minimum 

STRAG2 

DCCRCT 

CB 

100 

Branch current 



CL 

50 

Load current operating point 



CS 

12 

Source current operating point 



DELTA 


Tolerance on node voltage solutions 
(normally 10-5) 



ES 

12 

Source voltage operating point 



ITER 


Circuit solution internal iteration 
counter 



NITER 


Maximum allowable iterations in the 
circuit polution 



NOR 


Circuit reference node 



NSC 


Number of I-V points in Source (I) 
used in the circuit solution 



NT3S 


Number of points used in T3SRCS 



NT4S 


Number of points used in T4SRCS 



PP 


Constant power load for LOAD (I), 
variable not used when Phase II is 
driven by an interface tape 



PR 


Constant resistive power load for 
LOAD (I), PP(I) and/or PR(I) may be 
used to represent LOAD (I), variable 
not used when Phase II is driven by 
an interface load 



R 


Branch (I) line resistance 



RLOAD 

3 51 

Three point equivalent load 
I = load at 24v, 28v, 32 v 
J = 50 possible loads J = 51 voltage 
equivalent 





SEPS DATA VARIABLES 


PDP COMMON VARIABLE 

ELEMENT BLOCK NAME I_ J_ J<_ DEFINITION 

STRAG2 DCCRCT S TOO Branch (I) switch position 

(1 = closed, 0 = open) 

SC 15 2 12 Source curves 

l~ Number of points in curve 

J=1 Voltage 
J=2 Current 
K= Source number 



T3SRCS 

15 2 

Third source type 
1= Number of points in curve 

J=1 Voltage 
J=2 Current 


T4SRCS 

15 2 

Fourth source type 
1= Number of points in curve 

J=1 Voltage 
J-2 Current 


V 

30 

Node voltage 


V L 

50 

Load voltage operating point 


Z 

100 

Admittance of the branch 

CRTFLG 

IACSOL 

IDCSOL 


AC circuit solution flag 
DC circuit solution flag 

INVERT 

ACPOW 

9 

AC load for inverter AC BUS (I) 
variable not used when Phase II is 
driven by an interface tape 


INVOL 

10 

Inverter overload flag 


PFAC 

9 

AC load power factor for AC BUS (I) 


corresponding to loads in ACPOW(I). 
Not used with interface tape 


A- 8 



SEPS DATA VARIABLES 


PDP COMMON VARIABLE 

ELEMENT BLOCK NAME I_ J_ 

STRAG2 DEFCRT CRCT TOO 6 


ICHRBP 

I CRCT TOO 4 


INVLCA 9 

INVLCD 10 2 


LOADS 50 
NCRT 


K_ DEFINITION 

Branch (I) diode or RPC definition 

1= Branch No, 

J=1 Diode voltage drop 

<3=2 Diode forward resistance 

J=3 Reverse resistance 

J=4 RPC no load resistance 

J=5 Forward or reverse resistance 
being used 

<3=6 Branch current limit 

Load location of the battery charger 

Branch (I) definition 

1= Branch No. 

<3=1 Node number current out 

<3=2 Node number current in 

<3=3 Source number 

<3=4 Load number 

Inverter number connected to ac 
inverter bus (I) 

Table of branch no. vs inverter 
in the branch 

I* No. of entries in the table 

<3=1 Branch No. containing dc load 

J=2 Associated inverter no. for 

above dc load 

Branch locations of the loads 

Maximum or highest branch number 
in I CRCT 


A-9 



SEPS DATA VARIABLES 


POP 

ELEMENT 

STRAG2 


STRAG3 

STRAG4 


COMMON VARIABLE 


BLOCK 

NAME 

I_ J_ K_ 

DEFINITION 

DEFCRT 

NINVLA 


No. of inverters 


NINVLA 


No. of ac load buses supplied by 
i nverters 


NINVLD 


No. of dc loads used for inverters 


NLDS 


Highest load number 


NNODE 


Highest node number used in circuit 


NSCS 


Number of sources 


NSOUR 

12 3 

Source definition table 




1= Source number 




0=1 Branch location 


J=2 Source type 

0=3 Source reference numbers 

UV 30 Undervoltage limit for node (I) 

Not Used 


BATTRY AH 6 

BC 6 

BV 6 


CC 6 

CHRGLD 


CSUBD 6 

DQ 6 

EFF 6 


Actual number of ampere hours 
remaining 

Operating current point of the 
battery 

Operating voltage point of the 
battery 

Ampere hour capacity of the battery 

DC power required by battery charger 
when charging batteries 

Battery amp-hour capacity for 
battery (I) 

Amount of heat generated by the 
battery in watts 

Battery decimal efficiency during 
charge and discharge 
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SEPS DATA VARIABLES 


PDP 

ELEMENT 

COMMON 

BLOCK 

VARIABLE 

NAME 

I_ J_ 

K DEFINITION 

STRAG4 

BATTRY 

EFFAVR 


Average amp-hour efficiency during 
charging 



ICHRG 

6 

Charge flag to initiate charging 
of battery (I) 



IT 

5 2 

Battery temperatures used in SOCA, 





1= Temperature 





J= Battery group 



NSOCA 


Number of points used in SOCA 
for each I-V curve 



PI 


Constant used in the battery heat 
generation equation 



P2 . 


Constant used in the battery heat 
generation equation 



P3 


Constant used in the battery heat 
generation equation 



P4 


Constant used in the battery heat 
generation equation 



SOC 

6 

Battery (I) initial state-of-charge 



SOCA 

7 6 

2 Battery I-V curves versus temperature, 
and battery group 





1= Number of points 





J=1 Current 





0=2-6 Voltage at temperature of IT 





K= Battery group 1 or 2 (1 = group 

of 3 common batteries) (2 = 
group of common batteries) 

(per cell voltage) 



TB 

6 

Battery (I) temperature 



TD 

6 

Time of last major battery discharge 



TSS 


Battery steady state temperature 
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SEPS DATA VARIABLES 


PDP 

ELEMENT 

COMMON 

BLOCK 

VARIABLE 

NAME 

I_ J_ K_ 

DEFINITION 

STRAG4 

BATTRY 

XNC 

6 

Number of cells in battery I 

STRAG5 

FUSEL 

DAT 

5 

Time the fuel cell has been operating 
in decimal hours (MET) 



DTM 


Internal time step in FUCLTM 



FCCP 

5 

Fuel cell current 



FCHOL 


Fuel cell high temperature limit - 
heater turns "OFF" 



FCHTL 


Fuel cell redline limit - diagnostic 
warning 



FCLTL 


Fuel cell lower temperature limit - 
heaters turn "ON" 



FCT 

5 

Operating temperature of fuel cell (I) 



FCTA 

9 7 

Array containing fuel cell I-V 
curves as a function of temperature 





I s Number of points 





J=1 Current value (amps) 





J=2 Voltage level at each temperature 

thru of FCTN 

J=7 



FCTN 

6 

Temperatures associated with the I-V 
curves of FCTA 



FCWP1 

5 

Parasitic pump and logic loads for 
fuel cell (I) - constant power 



FCWP2 

5 

Parasitic heater cyclic load for 
fuel cell (I) constant resistance 



HPT 


Hydrogen purge time 



HR 


Hydrogen purge rate 



HUR 


Hydrogen use rate based on amp-hour 
requirements 



H2I 


Initial quantity of loaded in lbs 
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SEPS DATA VARIABLES 


PDP 

ELEMENT 

STRA65 


COMMON 

BLOCK 

VARIABLE 

NAME 


J K DEFINITION 

FUSEL 

H20T 


Total quantity of H?0 accumulated 
in lbs 


HDUM1 


Dummy filler variable 


HDUM2 


Dummy filler variable 


HDUM3 


Dummy filler variable 


H2T 


Initial quantity of hydrogen 
available 


KHTR 

5 

Fuel cell parasitic load flag 
Value 0 implies no parasitic load 
Value 1 implies parasitic load 


NFCTA 


Number of current points used in FCTA 


OPER 

5 

Flag indicating the on/off condition 
of the fuel cell 


OPT 


Oxygen purge time 


OR 


Oxygen purge rate 


OUR 


Oxygen use rate based on amp-hour 
requirements 


021 


Initial quantity of 0 % loaded in lbs 


02T 


Initial quantity of oxygen available 


PIH 


Interval between hydrogen purges 


PIO 


Interval between oxygen purges 


RES 

5 

Equivalent resistance of the fuel 
cell parasitic load at 28 VDC 


SSTVI 

10 

2 The T-I curve which the fuel cell 

follows as the fuel cell temperature 
reaches its steady state value 




0-1 * Temperature 




J=2 Current 


TMAXFC 


Maximum time step through feul cell 
thermal model 
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SEPS DATA VARIABLES 


POP COMMON VARIABLE 


ELEMENT 

BLOCK 

NAME 

I_ J_ K_ 

DEFINITION 

STRAG5 

FUSEL 

TPH 

5 

Time of last hydrogen purge for fuel 
cell (I) 



TPO 

5 

Time of last oxygen purge for fuel 
cell (I) 



WPR 


Fuel cell water production rate 

STRAG6 




Not Used 

STRAG7 




Not Used 

STRAGA 

PHAS1 

IADC 

250 2 

Activity dictionary 





1= Dictionary element 





J=1 Activity number 





J=2 Drum address 



IAN 


Number of dictionary entries 



ICDC 

750 2 

Component dictionary 


1= Dictionary element 
<J=1 Component number 
0=2 Drum address 


ICN 

Number of dictionary entries 

IPDC 750 

Procedure dictionary 


1= Dictionary element 


J=1 Procedure number 


J=2 Drum address 

IPN 

Number of dictionary entries 

MM 

Number of mission phases 

TABORT 

Mission elapsed time to end the 
simulation, default is 500 hours 
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SEPS DATA VARIABLES 


PDF 

ELEMENT 

COMMON 

BLOCK 

VARIABLE 

NAME 

I_ J_ K_ 

DEFINITION 

N/A 

ACPOWF 

PFEFF 


Inverter efficiency 



RESLOS 


Total system line loss 

N/A 

ALTERN 

I ACT 


Number of entries in the compact 
dictionary 



ICDCA 

500 3 

Compact component dictionary 





I- Number of entries 





0-1 Component number 





J=2 Corresponding drum address 





J=3 Relative location in the 
component dictionary 



ICDCB 

750 

Component usage count 

N/A 

BSLOC 

IBUSC 

500 

Component load assignment (COMPACT) 



ISYSC 

500 

Component system assignment (COMPACT) 

N/A 

BSLOCA 

IBUSC 

750 

Component load assignment 



ISYSC 

750 

Component systems assignment 

N/A 

CYCLIC 

I 


Number of entries 



MS 

100 

Cyclic mode 



NS 

100 

Cyclic number 



PERS 

100 

Cyclic period 



PONS 

100 

Cyclic percent on 



TS 

100 

Cyclic type 



ns 

100 

Cyclic time to start cycle 
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SEPS DATA VARIABLES 


POP 

ELEMENT 

COMMON 

BLOCK 

VARIABLE 

NAME 

I_ 

N/A 

DRMFLG 

ISF 


N/A 

FCYCL 

KF 


N/A 

INVEFF 

ACEFF 

3 

N/A 

INVPF 

AC P FAC 

9 

N/A 

MPSF 

ISF29 

50 



ISF36 

50 



ISF38 

50 



ISF40 

50 

N/A 

MTRANT 

TIMEM 

50 

N/A 

PRINT 

ICARD 




IPRNT 




I RESET 




ITAPE 

- 

N/A 

SUBSTM 

ILOC 

25 


ILOCN 

TITLE 6 25 

N/A TLINF IDA 

IDRM 
IEND 
IFIL 
IOUT 
IOUTM 
NWL 

N/A TOTPWR PWRTOT 


K_ DEFINITION 

Override control on illegal 
component turn off 

Type of entry indicator 

Inverter efficiencies 

AC load power factors 

Same as IFLG29 

Same as IFLG36 

Same as IFLG38 

Same as IFLG40 

Mission elapsed time at end of 
mission phase 

Card input read 

Input print required 

Simulation reset point 

Tape input read 

Subsystem number 

Number of subsystems 

Subsystem name 

Present drum address 

Drum full flag 

End of Phase I flag 

Output file number 

Output record number 

' Maximum number of output records 

Number of drum words left 

Total source power 
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